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Abstract 

This  research  produced  a  new  class  of  skin  friction  gages  that  measures  wall  shear 
even  in  shock  environments.  One  test  specimen  separately  measured  wall  shear  and 
variable -pressure  induced  moment.  Through  the  investigation  of  available  computational 
modeling  methods,  techniques  for  accurately  predicting  gage  physical  responses  were 
developed.  The  culmination  of  these  model  combinations  was  a  design  optimization 
procedure.  This  procedure  was  applied  to  three  disparate  test  conditions:  1)  short- 
duration,  high-enthalpy  testing,  2)  blow-down  testing,  and  3)  flight  testing.  The  resulting 
optimized  gage  designs  were  virtually  tested  against  each  set  of  nominal  load  conditions. 
The  finalized  designs  each  successfully  met  their  respective  test  condition  constraints 
while  maximizing  strain  output  due  to  wall  shear. 

These  gages  limit  sources  of  apparent  strain:  inertia,  temperature  gradient,  and 
uniform  pressure.  A  unique  use  of  bellows  provided  a  protective  shroud  for  surface 
strain  gages.  Oil  fill  provided  thermal  and  dynamic  damping  while  eliminating  uniform 
pressure  as  a  source  of  output  voltage.  Two  Wheatstone  bridge  configurations  were 
developed  to  minimize  temperature  effects  first  from  temperature  gradient  and  then  from 
spatially  varying  heat  flux  induced  gradient.  An  inertia  limiting  technique  was  developed 
that  parametrically  investigated  mass  and  center  of  gravity  impact  on  strain  output. 

Multiple  disciplinary  computational  simulations  of  thermal,  dynamic,  shear, 
moment,  inertia,  and  instrumentation  interaction  were  developed.  Examinations  of 
instrumentation  error,  settling  time,  filtering,  multiple  input  dynamic  response,  and  strain 
gage  placement  to  avoid  thermal  gradient  were  conducted.  Detailed  mechanical  drawings 
for  several  gages  were  produced  for  fabrication  and  future  testing. 
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1  Introduction/Motivation 


Why  develop  a  skin  friction  gage?  Aeronautical  engineers  since  the  time  of  Sir 
George  Cayley  have  sought  to  take  advantage  of  the  interaction  between  flow  and 
structure  to  permit  flight.  The  British  engineer  proposed  separate  aerodynamic  lift  and 
drag  forces  acting  on  an  aircraft.1  Over  the  past  100  years,  designers  have  successfully 
balanced  these  forces  against  thrust  and  weight  to  achieve  ever  increasing  performance 
growth.  Skin  friction,  also  called  wall  shear  (xw),  is  one  critical  piece  of  this  design 
puzzle.  The  accurate  measurement  of  this  parameter  can  literally  mean  the  difference 
between  a  design  that  can  and  cannot  fly.  Unfortunately,  the  complex  flows  often 
encountered  in  high-speed  vehicles,  such  as  in  the  isolator  section  of  a  scramjet,  produce 
shock  and  turbulent  effects  that  make  wall  shear  difficult  to  measure.  The  dramatic 
possibilities  of  such  complicated  flows  are  evidenced  in  Figure  1-1.  In  fact,  this  image  is 
from  a  test  conducted  to  determine  the  capabilities  of  the  new  class  of  gage  built  through 
this  research. 


Figure  1-1.  Shadowgraph  in  Virginia  Tech’s  supersonic  tunnel  with  Mach  4  nozzle. 

At  relatively  low  speeds,  the  only  means  of  imposing  flow  force  to  the  body  is 
through  two  surface  loads:  a  normal  force  and  a  tangential  force.  Respectively,  these 
forces  are  pressure  at  the  wall  (Pw)  and  wall  shear.  Even  at  supersonic  speeds  where 
wave  drag  and  shock  impingement  become  dominant,  the  flow  effects  are  still  translated 
to  the  structure  through  these  surface  loads.  At  hypersonic  speeds,  flow  ionization 
produces  a  magnetic  field  that  can  affect  ferrous  structures.  This  additional  force  is 
currently  being  researched  for  flow  control.  Finally,  heat  flux  can  produce  thermal 
loading  that  must  be  considered  in  thermal-mechanical  systems. 
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Of  these  loads,  what  impact  does  skin  friction  have  on  overall  system 
performance?  The  answer  depends  on  the  flow  regime  being  considered.  The  focus  of 
this  research  is  hypersonic  flight  and  specifically  measuring  skin  friction  in  supersonic- 
combustion  ramjet  or  scramjet  engines.  Although  the  skin  friction  coefficient  (Cf)  is 
typically  a  low  value,  around  0.002-0.005,  the  high  dynamic  pressure  (q)  experienced  in 
hypersonic  flight  makes  this  parameter  extremely  important.4  This  occurs  through  the 
relation  xw  =  Cf  q.  Additionally,  the  shock-rich  environment  can  make  measuring  wall 
shear  difficult.  Therefore,  let’s  examine  the  specific  impact  of  just  the  wall  shear 
component  of  drag  on  scramjet  performance. 

1. 1  Scramjet  Performance 

First  let’s  delve  into  the  scramjet  engine  to  see  how  it  works.  The  scramjet  is  a 
mechanically  simple  system.  A  sample  cut-away  view  is  shown  in  Figure  1-2.  This 
shows  a  cross-section  of  a  convergent-divergent  nozzle  with  a  combustor  in  the  middle. 
This  profile  can  be  rotated  about  an  axis  of  rotation  or  extended  into  the  page.  In  Figure 
1-2  separate  sections  of  the  engine  are  labeled.  We  will  follow  the  flow  path  from  left  to 
right.  First,  the  unavoidable  fore  body  shock  precedes  the  inlet.  The  flow  then  forms  an 
oblique  shock  train  in  the  isolator  section.  This  allows  a  desired  static  pressure  rise 
without  as  grave  a  total  pressure  loss  as  with  a  normal  shock.  This  maintains  a 
supersonic  flow  throughout  the  engine  allowing  for  higher  ram  pressure  than  a  ramjet. 
The  fuel  is  injected  in  the  isolator  section  for  mixing.  The  mixture  is  then  combusted  and 
expanded  through  the  aft-body  nozzle. 


Figure  1-2.  Scramjet  cross-section. 
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To  gain  a  better  appreciation  of  skin  friction  impact  on  scramjets,  it’s  useful  to 
examine  this  engine  type  against  other  engines.  Figure  1-3  shows  the  flight  regimes  of 
several  classic  engine  types  burning  molecular  hydrogen  as  fuel.  The  colored  diamonds 
portray  the  performance  of  a  typical  engine  for  each  class  using  respective  cycle  analysis 
for  verification.  The  engine  performance  characteristics  for  each  of  these  representative 
engines  are  detailed  in  Table  1-1. 5,6  Figure  1-3  clearly  illustrates  the  trade-off  between 
fuel  consumption  and  velocity.  Additionally,  the  higher  velocity  air-breathing  engines 
cannot  function  in  lower-speed  flight.  The  chemical  rocket  has  the  advantage  of  being 
able  to  fly  across  a  broad  spectrum  of  velocities,  but  at  the  expense  of  additional  oxidizer 
and  correspondingly  higher  take-off  weight. 
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Figure  1-3.  Flight  regimes  of  various  propulsion  systems  burning  H2. 

The  data  from  Table  1-1  provides  the  propulsion  engineer  with  an  apples-to- 
apples  comparison  of  system  performance.  These  values  are  the  same  as  quoted  in 
industry  by  engine  manufacturers.  Each  engine  color  in  the  table  corresponds  to  the  same 
diamond  color  in  the  plot.  Again  the  clear  trade-off  is  between  speed  and  fuel 
consumption.  If  you  want  to  go  fast,  you  have  to  bum  more  fuel.  Fortunately  with  a 
scramjet,  you  can  go  really  fast  without  burning  as  much  fuel  as  with  a  rocket  motor. 
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Table  1-1.  Engine  performance  parameters  for  3  air-breathing  engines  burning  H2. 


H2 

[Ramjet 

Scramjet 

F/mo(lbf-s/lb) 

72.354 

111.2 

63.35 

S  (lb/(lbf-s)) 

1.09E-04 

2.52E-04 

4.59E-04 

Isp  (s) 

|3971 

2177 

Vo  (km/s)  j 

jl .  1 5 

3.05 

r|o 

0.333 

0.495 

0.545 

r|th 

0.594 

0.762 

0.545 

r)P 

0.56 

0.65 

0.999 

The  performance  capabilities  of  the  scramjet  make  it  suitable  for  many 
applications.  One  such  use  led  to  design  work  on  NASA’s  single  stage  to  orbit 
Aerospace  Plane  shown  in  Figure  1-4.  Along  with  staging  for  low-earth  orbit,  the 
scramjet’s  high  velocities  make  it  a  strong  contender  for  long  range  missiles  that 
absolutely,  positively  must  destroy  a  target  within  seconds,  guaranteed.  Because  these 
engines  do  not  carry  an  oxidizer,  the  payload  can  be  greatly  increased  over  a  similarly 
classed  rocket.  Hopefully,  these  arguments  have  shown  that  scramjets  are  not  only 
intriguing  and  worthy  of  study,  they  actually  have  a  market  niche  that  is  currently 
untapped. 


:  -7— 

0 

•% 

Figure  1-4.  NASA’s  Aerospace  Plane  concept  art.7 


At  this  point  we  have  a  general  understanding  of  how  the  flow  moves  through  the 
scramjet.  In  order  to  understand  how  skin  friction  affects  performance,  we  must  apply 
some  rough  numbers  to  design.  The  dynamic  pressure  profile  through  the  engine  can  be 
approximated  using  the  stream  thrust  analysis  method.8  Figure  1-5  shows  this  dynamic 
pressure  variation  through  a  low-drag  engine.  We  will  examine  three  engines:  1)  Ideal  2) 
Low  Drag  and  3)  High  Drag.  Therefore,  Figure  1-5  shows  the  dynamic  pressure  profile 
through  Engine  2.  The  key  point  to  take  away  from  this  figure  is  the  high  dynamic 
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pressure  achieved  just  prior  to  the  combustor.  Even  with  a  low  skin  friction  coefficient, 
this  large  dynamic  pressure  produces  tremendous  drag. 


Figure  1-5.  Dynamic  pressure  profile  through  typical  scramjet  engine. 

Now  let’s  examine  the  performance  of  the  three  hypothetical  scramjet  engines. 
Engine  1  is  an  ideal  case,  where  skin  friction  is  assumed  to  be  negligible.  Engine  2  is  the 
low  drag  engine  of  our  dynamic  pressure  profile  where  the  wall  shear  in  the  isolator  and 
combustor  are  both  low,  but  reasonable  values.  Finally,  Engine  3  is  a  high  drag  case 
where  both  isolator  and  combustor  produce  large  wall  shear,  again  with  reasonable 
values.  A  real  engine  design  will  have  drag  characteristics  that  fall  between  the  Engine  2 
and  Engine  3  cases.  Figure  1-6  indicates  the  drag  coefficients  for  the  isolator  against  the 
combustor.  The  shaded  area  represents  a  bound  design  space  of  realistic  component  wall 
shear  values.9  This  broad  envelope  is  beyond  normal  expected  limits,  but  provides  a 
reasonable  starting  point. 

We  analyze  each  engine  under  the  same  flight  conditions  for  comparison. 

Realistic  free-stream  conditions  for  on-design  analysis  are  set  to  Moo  =10,  q«,=1000  lbf/ft2, 
and  Too=420  R.  Figure  1-6  shows  the  corresponding  T-s  diagram  for  each  engine.  The 
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relatively  high  entropies  reached  by  the  drag  engines  clearly  indicate  performance  losses. 
The  same  performance  parameters  used  to  compare  engines  of  different  type  delineate  the 
performance  of  each  of  our  scramjet  engines.  Table  1-2  provides  the  necessary  data  to 
determine  the  impact  of  wall  shear  on  engine  performance.  The  thrust  going  from  an 
ideal  to  a  low  drag  case  shows  a  6  %  reduction.  Even  more  severe,  the  high  drag  engine 
causes  a  31%  reduction  in  thrust  from  the  ideal  case. 


3  Engines 
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Figure  1-6.  Scramjet  engine  design  matrix  and  T-s  diagram:  1)  Ideal  2)  Low  drag  3)  High  drag. 


Table  1-2.  Performance  parameters  for  3  scramjet  engines. 


Engine  1 

Engine  2 

Engine  3 

F/mo  (lbf-s/lb) 

63.35 

59.65 

48.39 

S  (lb/(lbf-s)) 

4.59E-04 

4.88E-04 

6.01E-04 

Isp  (s) 

2176.99 

2049.81 

1662.85 

r|o 

0.545 

0.513 

0.416 

ijth 

0.545 

0.508 

0.398 

r)P 

0.999 

1.009 

1.046 

This  investigation  into  scramjet  performance  has  partially  answered  our  initial 
question.  Why  develop  a  skin  friction  gage?  Clearly  skin  friction  has  a  tremendous 
impact  on  overall  engine  performance.  We  have  seen  that  for  the  engine  type  we  are 
interested  in,  the  scramjet,  wall  shear  can  cause  between  a  6%  and  31%  reduction  in 
thrust  from  ideal.  The  designer  cannot  afford  to  guess  on  the  amount  of  drag  with  this 


6 


range  of  thrust.  Even  taking  the  conservative  route  by  assuming  a  3 1%  reduction  in 
thrust  requires  huge  structural  weight  increases  that  lead  to  other  problems.  But  why 
can’t  we  just  estimate  the  drag  from  computational  fluid  dynamics  (CFD)  or  some  other 
method?  This  is  the  topic  of  discussion  for  the  next  section. 

1.2  Estimating  Skin  Friction 

We’ve  seen  that  skin  friction  is  important,  but  can’t  we  just  estimate  it? 
Unfortunately,  a  CFD  turbulence  model  indirectly  relies  on  an  estimate  of  the  wall  shear. 
Other  techniques  exist  for  predicting  skin  friction  in  compressible  turbulent  flows.  The 
bad  news  is  that  for  each  technique  we  try  we  get  a  different  predicted  value.  Not  only  is 
this  frustrating,  but  the  only  way  to  resolve  the  issue  is  to  make  a  measurement.  With  a 
measured  value,  a  more  complicated  model  can  be  validated.  You  may  wonder  why  we 
even  make  a  model  if  we  have  to  make  a  measurement.  Actually,  one  only  need  measure 
the  wall  shear  at  a  point  as  an  anchor  for  a  complicated  model  where  the  wall  shear  varies 
across  the  surface. 

As  an  example  of  how  estimated  values  of  wall  shear  can  range  with  technique, 
let’s  examine  the  flow  field  in  the  Virginia  Tech  supersonic  tunnel  with  Mach  4  nozzle. 
Figure  1-7  shows  the  measured  boundary  layer  using  a  triple-rake  survey.  The  triple-rake 
is  located  0.046  m  downstream  of  the  Cf  gage  center  location.  This  boundary  layer 
information  is  used  directly  in  each  method  to  determine  skin  friction. 


Figure  1-7.  Boundary  layer  profile  Virginia  Tech  supersonic  tunnel. 
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Three  separate  techniques  of  estimating  wall  shear  are  employed:  Van  Driest  II, 
Schlichting,  and  Eckert. 10,1 1,12  The  Virginia  Tech  supersonic  tunnel  conditions  can  be 
modeled  as  steady  flat  plate  flow.  The  only  added  complications  come  from 
compressibility  effects  and  high  Reynolds  number,  i.e.  turbulence.  Both  of  these 
complications  are  considered  in  each  of  the  three  models  employed.  Table  1-3  shows  the 
different  estimates  for  wall  shear  for  the  Virginia  Tech  supersonic  tunnel. 


Table  1-3.  Estimates  for  the  Virginia  Tech  supersonic  tunnel. 


Van  Driest  II 

Schlicting 

Eckert 

Cf 

0.712E-3 

0.641E-3 

0.503E-3 

Xw  (Pa) 

74.2 

66.8 

52.4 

UV»!2) 

1856.3 

1671.5 

2198 

This  table  clearly  illustrates  the  problem  with  wall  shear  estimation.  Even  using 
the  exact  same  flow  conditions,  for  a  rather  simple  high-speed  flow,  the  techniques 
produce  results  that  at  the  extreme  vary  over  40%.  Additionally,  the  flow  measurements 
were  taken  approximately  0.046  m  behind  the  location  of  the  skin  friction  gage.  Even  in 
this  short  space,  boundary  layer  growth  would  suggest  that  the  wall  shear  drops  from  the 
skin  friction  gage  location  to  the  location  where  we  measure  boundary  layer  height. 

The  relatively  straight  forward  flow  in  the  supersonic  tunnel  is  easier  to  model 
than  in  real  scramjets  or  high  enthalpy  facilities  like  the  Large  Enthalpy  National  Shock 
(LENS)  facilities.  In  these  hot  reacting  flows,  simpler  transformation  models  break 
down  due  to  dissociation.  We  can  no  longer  apply  the  perfect  gas  model.  Fortunately, 
the  Eckert  reference  enthalpy  method  can  apply  to  these  flows.14  This  technique  is 
applied  to  the  conditions  expected  at  LENS  and  produces  a  reasonable  estimate  of  2000 
Pa  during  testing. 

The  preceding  arguments  have  been  based  on  two  assumptions  that  could  be 
completely  wrong  depending  on  the  location  of  interest  in  a  scramjet.  The  first  is  that 
there  is  a  skin  friction  to  measure.  This  assumption  could  be  completely  invalidated  in  an 
isolator  section  with  shock-impingement  that  has  caused  boundary  layer  separation.  In 
such  a  location,  the  wall  shear  might  be  zero,  or  the  flow  may  have  reversed  direction  in  a 
small  recirculation  zone  so  that  the  wall  shear  would  be  negative.  The  second 
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assumption  is  that  the  wall  shear  has  some  steady  value.  In  the  turbulent  flows  we’re 
examining,  the  flow  field  is  by  definition  unsteady,  and  this  translates  to  a  time  varying 
wall  shear.  As  followers  of  the  law  of  the  wall,  we  expect  that  any  variance  would  be 
minimal  within  the  laminar  sub  layer,  but  the  possibility  of  time  varying  wall  shear  exists. 
As  proof  of  the  relation  between  wall  shear  changes  and  turbulence,  consider  the  success 
of  Van  Driest’ s  oscillating  wall  model  in  representing  turbulence. 

Through  evaluation,  we  can  select  a  location  to  measure  the  wall  shear  that  is 
relatively  constant,  positive,  and  away  from  shocks.  Using  this  location  as  an  anchor  for 
a  complicated  CFD  model,  we  may  still  capture  much  of  the  complicated  flow  physics 
with  the  believability  that  in  at  least  the  relatively  calm  region  we  are  able  to  successfully 
match  a  measured  wall  shear  value. 

The  problem  arises  when  we  have  placed  a  gage  in  an  isolator  that  sees  off-design 
parameters.  This  is  the  reality  of  a  gage  used  in  flight  test.  Such  off-design 
characteristics  will  shift  the  shock  train  in  unavoidable  ways  that  might  manage  to  place 
an  oblique  shock  on  or  near  our  measurement  device.  If  we  mistakenly  ignored  such 
potential  problems,  we  could  report  wall  shear  values  that  are  completely  wrong.  This  is 
the  case  when  employing  a  robust  gage  incorrectly.  The  gage  designs  of  this  research 
provide  a  clear  method  of  avoiding  such  an  error. 

We’ve  answered  why  wall  shear  is  important  and  why  we  want  to  measure  it.  We 
need  to  decide  how  to  measure  it  and  how  we  know  whether  or  not  the  measurements  are 
reasonable.  This  is  the  subject  of  the  next  section. 

1.3  Measuring  Skin  Friction 

How  can  one  measure  just  the  wall  shear  in  a  flow  field  with  complex  turbulence, 
shock-impingement,  and  varying  heat-flux?  What  if  we  could  actually  separate  out 
pressure  variances  due  to  shocks  from  wall  shear  and  give  reasonable  measures  of  both 
the  shock  effect  and  the  wall  shear?  We’ve  jumped  a  little  bit  ahead,  but  this  is  exactly 
the  concept  that  is  proposed  and  proven  through  this  research.  Before  going  too  far, 
consider  methodologies  for  measuring  wall  shear. 
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There  are  a  variety  of  methods  for  measuring  skin  friction.  Each  of  these  has 
advantages  and  disadvantages  which  normally  have  to  do  with  the  flow  to  be  measured. 
Figure  1-8  provides  a  classification  hierarchy  for  methods  of  measuring  skin  friction. 


Figure  1-8.  Skin  Friction  measurement  technique  classification  hierarchy. 

Orr  provides  an  excellent  overview  of  many  methods  for  measuring  skin 
friction. 15  Of  these,  the  direct  measure  method  is  desired  for  these  high  speed  high  heat 
flux  environments.  A  tabulation  of  direct  measure  methods  to  date  is  available  at 
http://www.aoe.vt.edu/~ischetz/Skin  Friction/Cf  Gage  Table&Refs.pdf.  None  of  these 
previous  gages  provide  a  method  of  separating  wall  shear  and  pressure  variation  effect 
except  for  those  written  by  the  author.  This  is  a  new  area  of  research. 

The  initial  concept  for  this  new  kind  of  gage  has  its  basis  in  previously  built 
cantilever  beam  type  gages.  A  non-nulling  gage  is  preferred  for  simplicity  and  faster 
time  response.  Figure  1-9  shows  an  idealized  version  of  a  cantilever  beam  direct  measure 
gage  on  the  left.  In  this  gage  design,  the  flow  moves  across  the  head  surface  causing  a 
deflection.  This  deflection  can  either  be  measured  directly,  or  the  resulting  strain  can  be 
measured,  as  is  our  case.  One  advantage  to  this  design  is  its  ability  to  have  little  motion 
due  to  the  relatively  strong  wall  pressure  loading  while  producing  a  relatively  large 
output  due  to  the  small  wall  shear. 
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Figure  1-9.  Idealized  skin  friction  gage  cut-away. 


The  right  side  of  Figure  1-9  shows  a  3-D  rendering  of  the  bellows  type  gage 
developed  through  this  research.  It  contains  many  of  the  same  features  as  the  idealized 
skin  friction  gage.  The  main  similarity  is  the  exposed  surface  that  is  influenced  by  the 
flow  and  the  beam  that  flexes  to  the  load.  The  main  difference  is  that  the  bellows  forms  a 
cavity  that  traps  the  oil  filler  and  keeps  it  free  from  air  bubbles.  As  will  be  described  in 
this  paper,  the  oil  fill  provides  thermal  and  damping  characteristics,  but  more 
importantly,  it  helps  to  eliminate  uniform  wall  pressure  as  a  variable. 

1.4  Skin  Friction  Gage  Design 

A  general  assembly  cartoon  for  the  oil-filled  bellows  skin  friction  gage  is  shown 
in  Figure  1-10.  Here  we  see  the  main  components  that  make  up  this  gage.  The  first 
image  shows  the  element,  and  the  sensing  surface  is  highlighted.  The  head  and  root 
strain  gages  are  applied  at  the  locations  on  the  beam  as  shown.  The  stem  provides  a 
surface  for  solder  tabs,  strain  relief,  and  a  handle  for  proper  alignment.  In  the  next  image 
we  see  how  the  bellows  slides  up  and  around  the  element.  The  base  follows  suit,  and  all 
three  components  are  epoxied  together.  The  bellows  cavity  is  filled  with  oil,  and  a  plug 
is  inserted  to  prevent  leaks.  Finally,  this  entire  assemblage  is  placed  in  a  housing  with  a 
small  gap  at  the  head  for  motion. 
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Figure  1-10.  Assembly  of  oil-filled  bellows  skin  friction  gage. 


Through  this  research,  we  have  designed  four  separate  gages  for  fabrication,  CF1 
through  CF4.  These  gages  are  shown  in  Figure  1-11.  Each  gage  is  made  with  an  element 
which  is  green,  a  base  which  is  grey,  and  a  housing  which  is  pink.  Each  gage  also  has  a 
bellows  shown  in  blue.  The  critical  component  that  determines  how  the  gage  functions  is 
the  element.  This  component  has  the  exposed  sensing  surface  that  is  affected  by  the  flow 
and  the  beam  whose  strain  distribution  is  used  to  measure  wall  shear.  All  of  these  gages 
share  the  general  head  and  beam  of  the  idealized  gage  in  Figure  1-9,  even  though  the 
dimensions  are  different.  These  gages  vary  in  geometry  and  material.  CF2  through  CF4 
use  a  FC8  bellows  and  Lee-plug  for  oil  fill.  CF1  uses  a  small  FC3  bellows  without  oil 
fill.  CF1  also  lacks  a  stem.  All  gages  are  made  of  aluminum  except  CF3,  which  is  made 
of  INVAR. 
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CF1  CF2  CF3  CF4 


Figure  1-11.  Skin  friction  gages  from  this  research  to  scale. 

These  gages  are  design  solutions  to  measure  wall  shear  while  exposed  to  the 
nominal  set  of  load  conditions  shown  in  Table  1-4.  These  conditions  represent  a 
simplified  steady  state  for  a  high-enthalpy,  short-duration  test.  The  steady  state  for  a  test 
of  this  type  is  about  10  ms.  The  wall  shear  is  assumed  to  act  mainly  in  the  direction  of 
the  mean  flow  path  and  spanwise  shear  is  neglected.  The  acceleration,  ax,  is  caused  by 
the  motion  of  the  test  model  being  measured.  This  motion  is  also  dominant  in  the  mean 
flow  direction.  Finally,  the  heat  flux,  qw,  acts  perpendicular  to  the  sensing  surface. 

Table  1-4.  Nominal  design  specifications. 


Time 

(ms) 

T 

wx 

(Pa) 

ax 

(m/s2) 

(kW/m2) 

Load 

10 

2000 

9.8 

100 
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Each  gage  should  accurately  measure  wall  shear  and  be  insensitive  to  pressure, 
temperature,  and  shock  impingement.  Since  all  of  the  gages  of  this  research  were 
designed  to  these  nominal  specifications,  a  good  question  is  ‘why  are  they  so  different?’ 
These  gages  have  evolved  with  our  understanding  of  the  problem.  Though  each  gage  has 
successfully  met  our  requirements  at  each  step  along  the  way,  those  requirements  have 
been  finely  tuned.  The  culmination  of  these  efforts  is  a  design  optimization  methodology 
that  can  take  high  fidelity  requirements  and  provide  the  best  solution  to  meet  them.  In 
fact  CF3  and  CF4  designs  are  optimizations  to  the  highly  evolved  requirements 
developed  through  this  research. 

The  overall  objectives  of  this  research  follow: 

1 .  Design,  build,  and  test  a  gage  that  separately  measures  wall  shear  and  moment. 

2.  Determine  the  best  gage  characterization  techniques. 

3.  Develop  an  optimization  method  for  gage  design. 

4.  Use  the  optimization  technique  to  design  gages  for  three  disparate  test  conditions: 

1)  short-duration,  high-enthalpy  testing,  2)  blow-down  testing,  and  3)  flight  testing 

5.  Develop  simulations  for  complicated  physics. 

In  order  to  understand  the  design  methodology  of  the  bellows  direct  measure 
gage,  Chapter  2  presents  a  complete  idealized  simulation.  This  theoretical  approach  is  an 
overview  of  load,  gage,  and  instrumentation  interaction.  It  provides  the  foundation  for 
gage  geometry,  strain  gage  placement,  and  Wheatstone  bridge  design.  It  also  explains 
why  the  gages  shown  in  Figure  1-11  share  many  common  traits.  This  theoretical 
approach  is  directly  applied  to  achieve  the  first  objective. 

The  first  objective,  to  separately  measure  wall  shear  and  moment,  is  addressed 
through  the  design,  development,  and  testing  of  the  CF1  and  CF2  gages.  Chapter  3 
presents  the  evolutionary  process  and  experimental  results  for  CF1  and  CF2.  The  oil- 
filled  bellows  of  CF2  eliminates  uniform  pressure  as  a  variable  and  increases  damping. 
The  most  important  contribution  in  this  section  is  the  experimental  proof  of  the  ability  of 
CF2  to  separately  measure  wall  shear  and  moment  when  impinged  with  a  shock. 

Chapter  4  examines  the  performance  of  each  gage  to  the  nominal  load  conditions. 
Through  this  process,  objective  2  is  met.  The  best  method  for  accurately  predicting  gage 
response  is  determined.  We  select  the  best  model  to  predict:  1)  wall  shear  induced  strain 
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and  displacement  2)  shock  impingement  induced  strain  and  displacement  3)  vibrational 
frequency  4)  inertial  effects  and  5)  thermal  strain. 

Once  the  modeling  techniques  to  characterize  gages  have  been  built,  they  can  be 
used  to  develop  an  optimization  procedure  for  a  specific  set  of  load  conditions.  Chapter  5 
presents  the  development  of  this  procedure  meeting  the  3rd  objective.  It  also  provides 
optimized  solutions  for  the  three  test  conditions  listed  in  the  4th  objective.  This  is 
accomplished  through  the  development  of  the  CF3  and  CF4  gages.  First,  initial  decisions 
bound  the  geometric  design  space.  This  space  is  then  discretized  so  that  we  can  examine 
hundreds  of  possible  gages.  The  imposed  constraints  on  state  variables  further  collapse 
the  possible  solution  set.  Finally,  our  desired  performance  parameters  are  applied  to 
select  the  best  design  for  application. 

Chapter  6  uses  simulations  to  answer  the  tough  physical  questions  that  developed 
during  this  research.  These  higher  tiered  models  look  at  flow,  gage,  and  instrumentation 
interaction.  This  final  objective  allows  us  to  answer  the  following  questions:  What  filter 
should  be  applied?  How  will  the  gage  respond  to  a  time  varying  load?  How  about  time 
varying  inertia?  Then,  we  look  at  the  thermal  effects  on  measurement.  Could  errors  in 
strain  gage  placement  cause  a  voltage  output  due  to  temperature?  What  kind  of  error  can 
be  expected  from  instrumentation?  Could  a  spatially  varying  heat  flux  cause  apparent 
strain? 

Finally,  Chapter  7  presents  conclusions  and  recommendations.  This  chapter 
details  the  contributions  from  this  research.  It  also  discusses  how  well  we  met  our 
objectives.  Finally  it  provides  short  and  long  term  recommendations. 
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2  Theoretical  Approach 


The  gage  development  is  iterative,  but  must  start  with  a  theoretical  model.  A 
good  model  is  as  simple  as  possible,  captures  the  essence  of  the  physics,  and  produces 
useful  results.  In  our  case,  we  are  interested  in  accurately  modeling  the  flow,  the  gage 
response,  and  the  instrumentation  output.  Though  we  will  examine  complex  physical 
simulations  and  optimization  later,  this  chapter  presents  the  first  iteration  in  the  modeling 
process.  This  idealized  simulation  assumes  2-D  steady  flow,  a  cantilever  beam  gage,  and 
ideal  instrumentation.  This  low-fidelity  combination  provides  the  most  elegant  means  for 
understanding  the  design  features  of  this  new  class  of  skin  friction  gage.  It  unveils  the 
fundamental  means  for  separating  out  pressure  variations  from  wall  shear. 

2. 1  2-D  Steady  Flow  Model 

The  first  flow  model  for  simplicity  is  2-D  steady  flow.  Figure  2-1  shows  a  cutaway 
view  of  a  generic  supersonic  wind-tunnel  cross-section  as  well  as  the  flow-field 
coordinate  system  used  throughout  this  research.  The  blue  arrow  represents  the  mean 
flow  direction  which  is  parallel  to  the  x-axis,  and  the  y-axis  is  perpendicular  to  the  flat 
plate.  The  coordinate  system  is  centered  on  the  top  surface  of  the  skin  friction  gage. 


Floor  plate  ■  Beam 

H  □  oil 


Figure  2-1.  Cut-away  view  of  tunnel  test  section  with  flow-field  coordinate  system. 

In  this  2-D  steady  case,  only  the  mean  turbulent  flow  is  considered.  The  two 
surface  loads,  pressure  and  wall  shear,  act  on  the  structure.  Should  a  shock  impinge  on 
the  surface,  the  wall  pressure  distribution  varies  with  axial  length.  This  particular 
situation  will  be  explored  in  greater  detail  in  the  next  section.  The  heat  flux  and 
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acceleration  are  assumed  constant.  In  summary,  we  are  interested  in  four  simplified  load 
conditions  shown  in  Table  2-1.  Notice  that  wall  shear  and  inertial  loads  are  assumed 
axial. 

Table  2-1.  2-D  steady  flow  conditions. 


Wall  Shear 

Pressure 

Heat  Flux 

Acceleration 

Twx 

<*x 

2.2  Cantilever  Beam  Gage  Model 

The  previous  section  outlined  a  simple  set  of  flow  load  conditions.  The  simplest 
gage  model  available  is  the  cantilever  beam.  When  switching  to  the  gage  reference 
frame,  we  first  focus  on  the  axially  varying  pressure  distribution.  This  is  the  case  when  a 
2-D  shock  impinges  on  the  head  of  the  skin  friction  gage  as  shown  in  Figure  2-2.  This 
critical  case  is  encountered  in  many  high-speed  flow  situations  including  the  shock  train 
that  develops  in  scramjet  isolators.  Zooming  in  on  the  gage  itself,  we  can  begin  to 
examine  how  this  flow  condition  acts  on  the  gage. 


A  simple  gas  dynamics  analysis  of  the  flow  situation  presented  in  Figure  2-2 
breaks  the  flow  into  three  separate  regions.  The  flow  acting  on  the  front  half  of  the  gage 
surface  is  at  a  relatively  low  static  pressure,  Pi.  The  flow  acting  on  the  rear  half  of  the 
gage  surface  has  traveled  through  an  oblique  and  reflected  shock  and  is  at  a  higher  static 
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pressure  P3.16  This  2-D  distribution  is  shown  in  the  second  figure.  From  the  gage 
perspective,  this  pressure  distribution  has  an  equivalent  moment  and  pressure  load  as 
shown  in  the  third  figure.  These  are  M0  =  (P3-Pi)2rh3/3  and  P  =  (P3+Pi)Ah/2 
respectively.17  Since  P  is  simply  a  uniform  wall  load,  from  the  beam  reference  frame,  it 
can  be  applied  as  a  point  load.  We  will  rename  this  load  Pw  for  the  shock  impingement 
case. 

Now  we  can  describe  any  2-D  flow  load  acting  on  the  gage  head  as  three  simple 
separate  loads  from  the  beam  perspective.  The  beam  reference  load  conditions  are  shown 
in  Figure  2-3.  From  the  beam  reference  frame  a  complicated  axially  varying  pressure 
load,  Pw(x),  has  been  broken  into  two  simple  equivalent  loads,  Pw  and  M0.  The  wall 
shear  remains  the  same.  Notice  that  the  x  subscripts  have  been  dropped,  since  the  vectors 
define  the  direction. 


Figure  2-3.  Cantilever  beam  with  surface  flow  loads. 

A  convenient  effect  of  using  the  beam  reference  frame  is  that  each  of  these  loads 
is  well  understood  in  beam  theory.  For  example,  the  uniform  pressure  load  compresses 
the  beam  and  provides  a  constant  but  relatively  small  strain  along  its  length.  The  moment 
load  places  the  beam  front  face  in  tension  and  the  rear  in  compression.  The  wall  shear 
also  places  the  beam  front  in  tension  and  rear  in  compression.  In  the  case  of  the  moment 
load,  the  strain  produced  is  a  constant  for  each  side  of  the  beam.  However,  the  strain 
produced  by  wall  shear  is  a  maximum  at  the  beam  root  and  gradually  decreases  to  zero  at 
the  head.  The  effects  of  each  of  these  load  cases  are  shown  qualitatively  in  Table  2-2.  In 
this  table,  the  term  “Tension”  refers  to  the  upstream  face  of  the  beam  since  the  beam  is  in 
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tension  for  wall  shear  or  moment  loads  under  these  flow  conditions.  “Compression” 
describes  the  downstream  beam  face.  These  particular  images  were  developed  using 
finite  element  modeling  which  will  be  covered  extensively  during  gage  characterization. 

Table  2-2.  Strain  distributions  for  typical  beam  loads. 


Although  the  real  inertial  and  heat  flux  loads  are  quite  complicated,  we  can  take 
some  simplifications  from  our  simple  2-D  steady  flow  model.  Since  we  have  dictated 
that  the  inertial  load  is  purely  axial,  we  can  calculate  an  equivalent  acceleration  shear. 
Thus,  if  the  strain  produced  by  inertial  loading  is  much  smaller  than  that  from  wall  shear, 
we  can  ignore  inertia.  Since  the  heat  flux  is  assumed  constant,  any  temperature 
distribution  will  be  gradual  from  the  head  to  the  root  of  the  beam.  Therefore,  we  expect 
the  temperature  on  the  tension  side  of  the  beam  to  mirror  the  temperature  on  the 
compression  side.  This  is  a  useful  assumption  since  the  strain  varies  directly  with 
temperature. 
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The  fundamental  output  we  seek  for  any  gage  model  is  the  elastic  strain  along  the 
beam  axis  at  each  strain  gage  location.  In  other  words,  our  generalized  desired  output  is 
s(x,y,z,t).  We  are  interested  in  the  elastic  strain  in  line  with  beam  elongation  or 
shortening,  sy.  The  y  subscript  refers  to  the  flow  coordinate  system  y,  perpendicular  and 
up  from  the  floor  plate  and  aligned  with  the  strain  gages.  This  particular  strain  shows 
changes  along  the  beam  length.  This  strain  sums  the  outputs  of  five  sources  of  strain: 
wall  shear,  s™,  pressure,  epw,  thermal  gradient,  St,  moment,  smo,  and  inertial  loading,  ea. 
These  strains  correspond  to  the  loads  shown  in  Figure  2-4. 


Figure  2-4.  Load  configuration  for  simple  cantilever  beam. 


2.3  Ideal  Wheatstone  Bridge  Model 


The  gage  motion  produces  a  voltage  output  signal  through  strain  gages  placed 
along  the  beam  wired  through  Wheatstone  bridges.  As  each  strain  gage  elongates  or 
shortens,  its  resistance  changes  proportionally  through  Equation  2-1.  By  applying  a 
constant  excitation  voltage,  the  output  voltage  provides  a  direct  measure  of  the  strain. 
Strain  gages  allow  us  to  determine  beam  strains  at  a  variety  of  locations.  Wheatstone 
bridges  provide  the  opportunity  to  manipulate  these  outputs  in  a  useful  way. 


_  AL  _  1  AR 
£~  L~  GF  R 


Equation  2-1 18 
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A  configuration  using  two  half-bridges  consisting  of  two  strain  gages  per  bridge  is 
presented.  This  configuration  provides  output  results  from  measurements  at  the  root  and 
head  of  the  gage.  The  overall  ideal  Wheatstone  bridge  equation  is  presented  as  Equation 
2-2.  In  this  equation  for  each  of  the  half-bridges,  S3  and  84  are  zero. 

VGF  ■  Gain 


V.  =■ 


-U 


■  £2  +  £3 


,) 


Equation  2-2 
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2.4  Ideal  Simulation 

The  lowest  fidelity  model  allows  a  complete  examination  of  input  loading  to  gage 
response  to  expected  output  signal,  and  finally  to  verification  of  initial  load.  Figure  2-5 
shows  a  first  order  simulation  of  this  system.  The  load  conditions  acting  on  the  gage 
cause  changes  in  length  of  the  strain  gages  attached  to  the  surface.  The  heat  flux  causes  a 
temperature  variation  that  is  assumed  to  only  vary  down  the  axial  length  of  the  beam. 
With  this  model,  we  can  examine  the  expected  output  due  to  these  load  conditions. 


Figure  2-5.  First  order  simulation  of  load,  gage,  and  measurement  interaction. 

First  let’s  look  at  the  effect  of  a  uniform  wall  pressure  load.  Since  both  the 
tension  and  compression  side  strain  gages  see  the  same  compressive  strain,  we  know 
from  Equation  2-1  that  we  get  a  net  zero  voltage  output.  The  same  condition  occurs  due 
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to  a  temperature  distribution  that  is  only  a  function  of  vertical  direction.  Both  of  these 
load  conditions  are  imposed  on  the  system  in  Figure  2-6. 
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Figure  2-6.  Theoretical  output:  uniform  pressure  and  1-D  temperature  distribution. 

Since  we  have  eliminated  two  of  the  five  loads  from  Figure  2-5,  we  have  reduced 
our  loads  of  concern  to  three:  wall  shear,  moment  from  a  non-uniform  pressure,  and 
acceleration.  If  we  can  design  the  gage  to  minimize  acceleration  effects,  the  variables 
can  be  reduced  to  two:  wall  shear  and  moment.  This  theoretical  result  means  we  have 
simplified  our  problem  to  two  equations  and  two  unknowns.  Better  still,  these  loads  are 
linearly  independent,  so  we  have  reduced  our  problem  to  the  following  equations.  The 
four  constants  are  easily  found  through  calibration. 


Kr  =  CU^W+CUMo 

Equation  2-3 

VqH  =  C2l^w  C22 Mo 

Equation  2-4 

The  previous  developed  strain  gage  configuration  was  applied  to  CF 1  and  CF2 
and  is  going  to  be  used  for  several  of  the  CF3  and  CF4  gages.  Because  CF3  was 
designed  to  minimize  thermal-mechanical  effects,  we  do  not  expect  to  have  the  same 
thermal  problems  that  effected  CF2.  However,  there  is  a  theoretical  technique  to 
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eliminate  the  temperature  effects  even  if  the  temperature  varies  through  the  thickness  of 
the  beam.  This  situation  is  presented  graphically  in  Figure  2-7.  This  situation  might 
occur  if  the  heat  flux  varies  with  axial  flow  location. 


T 

T 


Figure  2-7.  Thermal  distribution  where  temperature  varies  through  beam  thickness. 


The  Wheatstone  bridge  configuration  applied  to  account  for  this  temperature 
variation  is  presented  in  Figure  2-8.  This  is  done  by  splitting  the  gage  along  its  centerline 
and  showing  the  tension  side  on  the  left  and  the  compression  side  on  the  right.  Each  side 
has  a  strain  gage  placed  horizontally.  This  horizontal  strain  gage  should  provide  little 
response  to  any  load  input  except  temperature.  Therefore,  the  thermal  effect  will  be 
cancelled  by  this  dummy  strain  gage.  This  strain  gage  configuration  will  be  tested  on 
some  of  the  CF3  and  CF4  gages. 
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Figure  2-8.  Two  full  Wheatstone  bridge  thermal  compensation  configuration. 
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3  Separately  Measure  xw  and  M0 

This  chapter  explains  the  design,  development,  and  testing  of  the  first  documented 
skin-friction  gage  to  separately  measure  wall  shear  and  moment.  This  chronology 
follows  the  development  of  the  CF1  gage  and  the  evolutionary  changes  that  made  CF2 
successful.  The  CF 1  gage  provided  the  first  work  using  the  bellows  as  a  thermal  barrier. 
The  CF2  gage  proved  the  concept  of  separately  measuring  wall  shear  and  shock-induced 
moment. 

3. 1  CF1  Gage  Design 

The  number  of  possible  design  choices  and  trade-offs  for  the  first  iteration  seemed 
overwhelming.  Fortunately,  thanks  to  the  theoretical  development,  we  neglected  certain 
state  variables,  albeit  somewhat  naively.  Table  3-1  summarizes  the  effects  considered 
and  approach  taken.  From  theory,  temperature  and  uniform  pressure  effects  were 
eliminated.  Next,  making  the  gage  as  light  as  possible  minimized  inertial  loading. 
Intuitively,  we  would  like  to  minimize  the  gap  size.  However,  under  low  wall  shear 
conditions,  the  gage  motion  was  determined  to  be  relatively  small.  Thus  displacement, 
hence  gap  size,  was  left  to  be  checked  after  design.  Moment  effects  were  also  left  for 
post  design.  This  limited  the  design  state  variables  to  wall  shear  strain,  frequency,  and 
damping.  Structural  damping  was  assumed  to  be  about  0.02  and  possible  methods  for 
increasing  this  value  were  considered  for  CF1.  These  included  using  grease  on  the 
bellows  and  tape  on  the  beam  for  viscous  damping.  Without  knowing  the  total  system 
impact,  damping  was  temporarily  shelved  as  a  design  variable.  Therefore,  the  goal  was 
to  achieve  the  most  strain  output  and  highest  frequency  possible.  Unfortunately,  these 
two  characteristics  oppose  each  other. 

Even  though  strain  and  frequency  have  become  the  focus  of  optimization, 
temperature  and  acceleration  considerations  really  drove  CF1  design.  In  order  to  protect 
the  strain  gages  from  convective  heat  transfer  while  allowing  nearly  unconstrained 
motion  of  the  overall  system,  we  wanted  to  use  a  bellows.  Working  with  commercially 
available  bellows  we  wanted  the  smallest  bellows  available  to  minimize  weight.  Thus  we 
chose  Servometer  ®  FC-3.  This  decision  greatly  constrained  the  amorphous  design  by 
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limiting  the  geometric  design  space.  Immediately,  the  head  thickness  and  head  radius 
were  set.  Additionally,  the  maximum  beam  radial  dimension  was  constrained  and  the 
beam  length  was  set.  It  also  made  optimization  of  strain  and  frequency  manageable. 
Thus,  our  design  variables  have  been  limited  to  material  choice  and  beam  cross-section. 


Table  3-1.  CF1  design  approach. 


Physics 

State  Variables 

Approach 

Optimize 

Temperature 

£at 

Neglect  (Theory) 

N/A 

Pressure 

Spw 

Neglect  (Theory) 

N/A 

Acceleration 

£a 

Lightweight 

Minimize  m(mass) 

Displacement 

8 

Post-Design  Check 

N/A 

Moment 

SMo 

Post-Design  Check 

N/A 

Wall  Shear 

£tw 

FEM 

Maximize  sxw 

Frequency 

f 

FEM 

Maximize  f 

In  order  to  maintain  a  reasonable  design  space,  the  materials  were  limited  to 
aluminum,  copper,  and  polyethersulfone,  a  high-temperature  plastic.  These  materials  had 
been  used  in  building  previous  gages.  Copper  is  a  typical  skin  for  scramjets,  aluminum  is 
easy  to  manufacture,  and  polyethersulfone  shows  little  thermal  response.  The  beam 
cross-sections  were  limited  to  circular,  rectangular,  and  annular.  Each  potential  design 
was  built  as  a  3-D  Baseline  FEM.  Figure  3-1  shows  the  resulting  response  data  as  a 
function  of  material  and  cross-section. 
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Material 


*  Circle  max  tensile  strain 

■  Rectangle  max  tensile  strain 
o  Annulus  max  tensile  strain 

#  Circle  Frequency 

■  Rectangle  Frequency 
o  Annulus  Frequency 

Material  1-  Copper 
Material  2-  Aluminum 
Material  3-  PES 


Figure  3-1.  Strain  and  frequency  as  a  function  of  material  and  cross-section. 
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This  figure  emphasizes  the  trade-off  between  vibration  frequency  and  strain.  The 
annular  design  posed  too  many  manufacturing  difficulties  to  achieve  an  exact  mass 
distribution.  In  other  words,  drilling  a  hole  exactly  at  the  center  is  difficult.  Also,  in 
order  to  apply  the  strain  gages  on  a  flat  surface,  the  circular  cross-section  was  discarded 
in  favor  of  the  rectangular  one.  After  selecting  a  workable  cross-section,  we  needed  to 
choose  the  best  material.  Because  the  structural  damping  coefficient  was  expected  to  be 
low,  we  decided  that  aluminum  provided  a  reasonable  balance  between  vibration 
frequency  and  strain. 

Once  an  aluminum  rectangular  design  was  selected,  the  exact  cross-section  area 
needed  specifying.  This  was  partially  limited  by  the  strain  gage  selection.  The  strain 
gages  were  selected  to  have  the  minimum  design  length,  aluminum  self-temperature 
compensation,  and  a  minimum  350  Q  gage  resistance,  Rg.  A  commercially  available 
Vishay®  EA- 13-03 ICE-350,  met  the  requirements  with  a  gage  length  of  0.076  in.  and 
width  of  0.062  in.  This  set  a  minimum  beam  width  for  gage  placement.  Now  the  beam 
width  was  constrained  to  a  strain  gage  minimum  and  bellows  maximum.  Then  a  sweep 
was  made  of  design  cross-section  effect  on  frequency  output.  Figure  3-2  shows  the  non¬ 
linear  relationship.  The  thinnest  tested  design  provided  a  reasonably  high  vibration 
frequency.  This  became  the  design  configuration  that  was  built  for  CF1. 


Figure  3-2.  CF1  frequency  sweep  as  a  function  of  beam  cross  section. 
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3.2  CF1  Results 


Each  CF1  Wheatstone  half-bridge  was  completed  with  a  Vishay®  Model  2310 
signal  conditioner.  These  amplifiers  were  set  using  the  10  Hz  low-pass  Butterworth 
filter,  a  5  V  input  voltage,  and  10k  gain.  A  National  Instruments  ®  data  acquisition 
system  was  controlled  by  programs  written  in  Lab  VIEW  ®.  The  sample  rate  was  500 
Hz.  Prior  to  each  test  set,  both  23 10’s  were  zeroed. 

The  CF 1  gage  met  many  of  our  design  criteria  requirements  during  static 
calibration.  Both  the  static  wall  shear  and  static  moment  calibration  had  high  linearity 
between  load  and  voltage  output.  CF1  also  exhibited  the  expected  difference  between 
wall  shear  induced  and  moment  induced  strain  distribution.  Table  3-2  shows  pictures  of 
the  gage  in  several  of  its  assembly  stages.  In  (a)  the  element  with  attached  strain  gages  is 
shown.  Next  (b)  shows  the  wires  attached.  The  picture  in  (c)  shows  the  element 
shrouded  with  bellows  and  base  attached.  Finally  (d)  shows  a  scaled  drawing  of  the 
assembled  housing,  element,  and  base. 

Table  3-2.  CF1  gage  in  stages. 


3.2.1  CF1  Calibration 

The  CF1  gage  underwent  static  wall  shear  and  moment  calibrations.  A  shaker- 
test  to  determine  vibration  frequency  was  attempted,  but  broke  one  of  the  two  CF1  gages. 
After  tunnel  run  testing,  thermal,  pressure,  and  wire  stiffness  sensitivity  tests  were 
conducted.  These  will  be  presented  after  experimental  test  results.  Once  the  pressure 
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issue  was  uncovered,  attempts  at  further  CF1  frequency  characterization  were  discarded 
in  favor  of  redesign. 


CF1  Static  Wall  Shear  Calibration 

The  static  wall  shear  calibration  was  conducted  by  placing  the  CF 1  beam 
perpendicular  to  earth’s  gravity  and  applying  weights  to  simulate  wall  shear.  This 
process  is  sketched  in  Figure  3-3.  A  string  was  taped  to  the  top  surface  with  an  attached 
low  weight  “bucket”.  The  string,  bucket,  and  tape  provided  minimal  weight,  but  were 
accounted  for  through  zeroing  of  the  signal  conditioners  prior  to  adding  weights. 


riamn 


Tape 

/-String 

/Weight 


Figure  3-3.  Static  wall  shear  calibration  technique. 

The  completely  fabricated  CF1  gage  was  calibrated  using  this  technique.  Several 
tests  were  conducted  varying  the  weight  and  order  of  weight  placement.  No  hysteresis 
effect  was  detected.  The  resulting  relationship  between  wall  shear  and  output  voltage 
was  highly  linear  as  evidenced  in  Figure  3-4. 
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Figure  3-4.  CF1  static  wall  shear  calibration  results. 

CFl  Static  Moment  Calibration 

The  static  moment  calibration  was  conducted  by  placing  the  CFl  beam  parallel 
with  gravity,  taping  a  lever  arm  to  the  gage  top  surface,  and  applying  weights  to  produce 
moment.  This  technique  is  sketched  in  Figure  3-5.  The  lever  arm  was  a  long  thin 
cylindrical  bar  that  was  centered  on  the  head  to  produce  zero  moment  without  weight. 
Additionally,  the  bar  was  placed  in  line  with  axial  flow  to  produce  a  direct  moment. 
Weights  were  then  added  at  a  constant  distance  from  the  head  center. 


Figure  3-5.  Static  moment  calibration  technique. 
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A  fully  fabricated  CF1  gage  was  calibrated  using  this  technique.  As  with  the  wall 
shear  calibration,  no  hysteresis  effect  was  noticed  despite  varying  weight  loading  and 
unloading  order.  The  high  linearity  between  loading  and  voltage  output  is  shown  in 
Figure  3-6. 


Moment  Calibration 


y=904.28x  +  0.0002 
R2  =  0.997 


♦  Upper(V) 

■  Root  (V) 

- Linear  (Upper(V)) 

- Linear  (Root  (V)) 

y=657.43x  + 0.0006 
R2  =  0.9947 


Figure  3-6.  CF1  static  moment  calibration. 


3.2.2  CF1  Virginia  Tech  Supersonic  Tunnel  Wall  Shear  Test 

The  validation  experiment  for  CF1  was  conducted  in  the  Virginia  Tech  supersonic 
test  facility  with  the  Mach  4  nozzle  installed.  This  is  the  blow-down  facility  shown  in 
Figure  3-7.  The  test  section  has  windows  located  on  either  side  of  the  tunnel  for  taking 
shadowgraphs.  Accelerometer  measurements  confirmed  that  this  tunnel  experiences 
negligible  acceleration. 

Unfortunately,  the  CF 1  experimental  test  results  could  not  be  correlated  to  wall 
shear,  static  temperature,  or  static  pressure.  With  so  many  possible  sources  of  apparent 
strain  and  the  possibility  of  instrumentation  error,  an  investigation  of  CF1  behavior  was 
conducted  post  experiment. 
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Figure  3-7.  Virginia  Tech  supersonic  tunnel  diagram.22 

3.2.3  CF1  Post  Experiment  Investigation 

During  validation  testing,  the  gage  exhibited  voltage  output  error  that  led  to  an 
investigation  of  sources  of  apparent  strain.  Each  of  these  possible  sources,  thermal, 
pressure,  and  wire-stiffening,  were  investigated.  All  of  these  provided  design  inputs  for 
CF2.  The  pressure  issue  was  dealt  with  by  using  incompressible  oil  in  the  bellows.  The 
manufacturing  concerns  were  the  need  for  a  hold,  strain  relief,  and  solder  pad  surface. 

All  of  these  were  solved  by  adding  a  handle  and  securing  the  wires  to  the  handle.  The 
thermal  issue  was  approached  by  using  the  oil  as  an  insulator  for  the  strain  gages. 

Finally,  wire-stiffening  was  solved  by  moving  the  wire  holes  in-line  rather  than 
perpendicular  to  the  beam. 

After  the  first  set  of  test  runs  in  the  Virginia  Tech  supersonic  tunnel  facility,  it 
was  clear  that  some  source  of  apparent  strain  was  affecting  the  voltage  output  for  CF 1 . 

In  order  to  either  account  for  this  error,  or  eliminate  it  through  design,  several  possible 
physical  sources  were  examined.  Static  pressure,  static  temperature,  and  wire  stiffening 
all  proved  to  be  possible  sources  of  apparent  strain. 

CF1  Vacuum  Test 

The  possibility  that  static  pressure  was  a  source  of  apparent  strain  required 
investigation.  Figure  3-8  shows  the  vacuum  test  setup.  The  CF1  gage  was  connected  to  a 
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vacuum  pump  with  a  vacuum  hose  secured  to  the  housing.  A  tee  was  used  to  connect  to 
a  50  psi  pressure  transducer  to  measure  static  pressure.  Additionally,  a  Type-T 
thermocouple  was  placed  near  the  head  of  the  CF1  gage  to  measure  temperature.  The 
other  leg  of  the  tee  connected  to  a  valve  and  vacuum  pump. 


Figure  3-8.  Vacuum  test  technique. 


This  test  was  conducted  with  CF1  at  multiple  angles  to  gravity.  The  vacuum  was 
applied  and  valve  closed  to  check  for  leaks.  The  gage  reaction  at  the  head  and  the  root  is 
shown  in  Figure  3-9.  Here,  the  white  trace  shows  the  changing  static  pressure 
measurement.  Both  the  head  output  voltage  in  red  and  the  root  output  voltage  in  green 
trend  very  closely  with  pressure.  The  temperature  changed  less  than  1  deg  C  during  this 
test. 

The  observed  pressure  effect  on  CF 1  measurements  was  intolerable.  A  possible 
solution  to  this  source  of  apparent  strain  was  to  fill  the  bellows  cavity  with 
incompressible  fluid.  By  doing  so,  the  volume  inside  the  bellows  would  remain  constant. 
This  would  negate  the  beam’s  ability  to  expand  or  contract  due  to  changes  in  uniform 
static  pressure.  This  design  solution  was  applied  to  the  CF2  gage. 
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CF1  Thermal  Test 

The  changing  static  temperature  measured  during  the  experiment  required 
examination.  A  temperature  test  was  conducted  with  the  CF1  and  a  thermocouple  placed 
in  a  small  oven  where  temperature  conditions  were  well  controlled.  Both  instruments 
were  packed  in  ice  and  allowed  to  warm  throughout  the  full  test.  Both  the  head  and  root 
voltages  varied  inversely  with  increasing  temperature  as  shown  in  Figure  3-10.  This  is 
shown  by  the  decreasing  red  and  green  traces  represented  by  CF1  voltages.  The 
thermocouple  output  is  shown  in  blue.  The  high- lighted  region  of  interest  represents  the 
temperature  range  encountered  to  the  tunnel  experiments. 

This  thermal  test  proved  temperature  was  a  source  of  apparent  strain.  The 
solution  of  the  pressure  problem  by  oil-fill  posed  a  possible  solution  to  the  thermal  issue. 
By  filling  the  bellows  cavity  with  oil,  it  was  hoped  that  the  temperature  effect  could  be 
significantly  reduced.  This  was  the  approach  taken  in  developing  CF2. 
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Figure  3-10.  CF1  post  experiment  thermal  test. 

CF1  Wire  Stiffening  FEM 

To  illustrate  the  impact  of  wire  placement  on  the  overall  strain  output,  examine 
Figure  3-11.  The  tension  side  shows  high  strain  at  the  wire  attachment  points  indicated 
by  the  gray  and  red  contours.  On  the  compression  side  figure,  high  compression  strain 
develops  at  the  wire  attachment  points  indicated  by  the  blue  contours.  Though  the  wires 
were  small  diameter,  they  were  also  short  and  stiff.  The  wire  stiffness  reduced  the 
amount  of  strain  that  would  normally  be  produced  by  xw.  Calibration  can  reduce  this  as 
an  issue  in  a  small  wall  shear  regime,  where  the  wires  act  elastically.  However,  we  want 
as  much  strain  produced  as  possible  at  the  strain  gage  location.  Therefore,  by  moving  the 
wires  to  the  sides  of  the  beam,  they  do  not  counteract  the  wall  shear. 
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Figure  3-11.  CF1  wire  truss  stiffening. 


3.3  CF2  Gage  Design 

The  CF2  design  scheme  was  nearly  identical  to  that  of  CF 1 .  Several  design 
changes  were  implemented.  The  focus  was  to  eliminate  strain  due  to  uniform  wall 
pressure.  Application  of  an  incompressible  fluid  within  the  bellows  would  not  permit  the 
beam  to  contract  or  elongate  under  uniform  pressure  load.  Additionally,  the  oil  benefit  in 
diffusing  thermal  flux  and  increased  damping  were  seen  as  side  benefits.  Therefore,  oil- 
fill  became  the  priority.  Because  of  the  geometry  required  to  accommodate  the  needed 
Lee  Plug®,  a  larger  bellows,  Servometer  ®  FC-8,  was  chosen.  This  was  the  smallest 
bellows  available  for  the  oil  configuration.  Keeping  the  gage  small  was  aimed  to 
minimize  inertial  effects.  As  an  aside,  another  technique  for  smaller  bellows 
implementation  with  oil  fill  is  proposed  in  the  recommendations  chapter.  Returning  to 
this  discussion,  the  FC-8  bellows  froze  the  head  geometry.  From  the  arguments 
presented  for  CF1,  we  stuck  with  the  rectangular  aluminum  cross-section. 

By  using  the  same  strain  gages  as  the  previous  design,  we  removed  Rg  as  a 
variable.  As  with  the  previous  optimization,  gap  size,  8,  and  acceleration  output,  sa,  and 
moment  output,  sm0,  were  left  for  verification  after  optimization.  Therefore  we  have 
narrowed  our  optimization  variables  to  those  in  Table  3-3.  The  biggest  difference 
between  the  CF1  optimization  and  this  was  the  realization  that  the  beam  length  could  be 
varied  below  the  limit  of  the  bellows.  This  additional  design  variable  meant  we  could 
vary  the  beam  length,  L,  width,  w,  and  thickness,  th. 
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Table  3-3.  CF2  state  variables,  objective  functions,  and  design  variables. 


State  Variables 

Maximize 

Minimize 

Design 

Strain  due  to  Wall 

8S 

N/A 

Ac,  L 

Frequency 

©n 

N/A 

Ac,  L 

To  limit  the  beam  length  and  thickness,  we  recognized  manufacturing  limitations. 
After  discussions  with  the  shop,  the  beam  length  was  set  to  0.475  in.  and  the  minimum 
thickness  for  the  beam  was  set  to  0.040  in.  With  these  geometric  constraints  in  place,  we 
could  vary  the  thickness  to  find  the  most  strain  output  that  still  provided  a  minimal 
natural  frequency.  This  turned  out  to  be  the  minimum  allowable  thickness. 

3.4  CF2  Results 

The  CF1  apparent  strain  investigation  provided  several  design  improvements  to 
CF2.  Once  these  design  characteristics  were  applied,  the  gage  was  built,  calibrated,  and 
had  separate  test  verifications  in  the  Virginia  Tech  Supersonic  Tunnel.  Table  3-4  shows 
several  stages  of  the  CF2  assembly:  a)  element  with  strain  gages  b)  3-D  rendered 
assembly  and  c)  scaled  drawing  comparison  against  CF 1 . 

Table  3-4.  CF2  gage  in  stages. 


Each  CF2  Wheatstone  half-bridge  was  completed  with  a  Vishay®  Model  2310 
signal  conditioner.  These  amplifiers  were  set  using  the  10  Hz  low-pass  Butterworth 
filter,  a  5  V  input  voltage,  and  a  10k  gain.  A  National  Instruments  ®  data  acquisition 
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system  was  controlled  by  programs  written  in  Lab  VIEW  ®.  The  sample  rate  was  500 
Hz.  Prior  to  each  test  set,  both  23 10’s  were  zeroed. 


3.4.1  CF2  Calibration 

The  gage  characterization  required  several  calibration  tests.  After  the  testing 
calibration  results  from  CF1,  a  calibration  matrix  was  developed  to  determine  sensitivity 
of  gage  performance  to  configuration.  The  three  tested  gage  configurations  were 
designated  Baseline,  Bellows,  and  Oil.  Respectively,  these  configurations  consist  of  the 
unshrouded  element,  element  with  base  and  bellows  attached,  and  finally  the  oil-filled 
complete  gage.  The  reason  for  calibrating  the  gage  under  these  different  configurations 
was  to  determine  the  effect  of  the  bellows  and  oil  on  gage  performance.  As  expected, 
there  was  little  or  no  effect  on  static  calibrations  through  the  different  configurations,  but 
a  significant  effect  during  frequency  characterization. 

CF2  Static  Wall  Shear  Calibration 

The  static  wall  shear  calibration  of  CF2  was  conducted  exactly  as  with  CF 1 .  An 
example  of  such  testing  is  shown  in  Figure  3-12.  This  particular  photograph  shows  the 
calibration  as  conducted  for  the  Baseline  configuration.  However,  the  results  for  the  Oil 
or  assembled  configuration  are  presented  in  Figure  3-13.  The  results  for  both  the  head 
and  root  voltage  output  show  high  linearity  with  wall  shear. 


Figure  3-12.  CF2  baseline  static  wall  shear  calibration  setup. 
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Figure  3-13.  CF2  static  wall  shear  calibration  results. 


CF2  Static  Moment  Calibration 

The  static  moment  calibration  was  conducted  as  presented  for  CF1.  Figure  3-14 
shows  the  set-up  for  the  Baseline  configuration  calibration.  The  calibration  curves  for  the 
head  and  root  shown  in  Figure  3-15,  meet  expectations.  Both  curves  are  highly  linear. 
Additionally,  the  proportionality  constant  for  the  head  and  root  are  nearly  the  same, 
within  7%,  as  expected  with  the  constant  strain  distribution  expected  for  a  pure  moment 
load. 


Figure  3-14.  CF2  baseline  static  moment  calibration  setup. 
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Static  Moment  Calibration 


▲  Root  + 

A  Head  + 

□  Root  - 

□  Head  - 

- Linear  (Root-) 

- Linear  (Head  -) 


Figure  3-15.  CF2  static  moment  calibration  results. 


CF2  Vacuum  Test  Calibration 

The  CF2  vacuum  test  was  the  most  important  calibration  to  ensure  that  our 
improved  design  eliminated  the  pressure  problem  of  CF 1 .  This  particular  calibration  was 
only  conducted  on  the  assembled  Oil  configuration.  Figure  3-16  shows  the  result  of  this 
calibration  test.  The  white  trace  is  the  pressure  transducer  output  throughout  the  test. 

The  red  and  green  traces  are  CF2  root  and  head  output  during  the  test.  The  test  order  was 
set  with  step  4  allowing  us  to  check  for  leaks.  Despite  the  small  amount  of  noise 
observed  throughout  the  test,  neither  set  of  voltage  output  shows  signs  of  pressure 
reaction.  This  is  especially  obvious  comparing  the  CF2  results  to  those  of  CF1. 
Therefore,  this  calibration  test  has  successfully  shown  that  a  uniformly  distributed 
pressure  on  the  top  surface  of  CF2  will  not  affect  voltage  output.  This  satisfied  our  most 
important  design  criteria  in  moving  from  CF1  to  CF2. 
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Figure  3-16.  CF2  vacuum  test  calibration. 


CF2  Thermal  calibration 

Several  tests  for  thermal  calibration  were  conducted.  Cold  and  hot  soaks  as  well 
as  gradually  varying  temperature  tests  were  conducted.  Though  all  tests  showed  at  least 
some  transient  periods  of  apparent  thermal  strain,  no  direct  relationship  between 
temperature  and  voltage  output  could  be  developed. 

CF2  Dynamic  response  calibration 

The  dynamic  response  information  we  desire  is  the  settling  time.  The  method  for 
test  was  a  step-release.  This  test  was  conducted  by  tying  a  thread  to  the  head  of  the  gage, 
imparting  a  simulated  axial  wall  shear,  cutting  the  string,  and  measuring  the  natural 
response.  A  sketch  of  this  process  is  shown  in  Figure  3-17.  The  CF2  gage  was  calibrated 
under  three  different  configurations.  The  Baseline  configuration  is  the  simply  the 
unshrouded  element.  The  Bellows  configuration  has  the  bellows  placed  on  the  element, 
but  not  glued  and  without  oil-fill.  Finally,  the  Oil  configuration  is  the  full  test 
configuration  with  oil  fill. 
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Figure  3-17.  Dynamic  response  calibration  setup. 


Figure  3-18  shows  the  results  of  the  string  tests.  These  results  are  presented  for 
quick  reference  in  Table  3-5.  Three  configurations  were  measured  for  parametric 
analysis  and  FEM  validation.  First,  from  a  parametric  point  of  view,  each  case 
completely  changes  the  dynamic  response.  Perhaps  the  most  important  characteristic  for 
short  duration  testing  is  the  settling  time.  For  example,  the  CF2  Oil  configuration  will 
not  settle  in  time  to  accurately  measure  xw  during  a  5  ms  test  run.  The  response  of  the 
bellows  case  should  not  be  weighed  too  heavily  since  the  components  were  not  glued 
together,  so  chatter  was  an  issue.  The  oil  configuration  also  showed  signs  of  second 
mode  excitation.  However,  this  does  not  present  an  issue  in  terms  of  determining  settling 
time. 


Figure  3-18.  String  test  calibration  data.  Top  row  shows  calculated  oo„  and  C.  Bottom  row  shows 
measured  rise  and  settling  times  for  +10%  load. 
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Along  with  the  presented  rise  and  settling  times,  an  estimation  of  vibrational 
frequency  and  damping  coefficient  for  each  configuration  is  given.  These  are  not 
necessary  in  characterizing  the  CF2  dynamic  response.  They  are  critical,  however,  in 
providing  a  means  for  future  modeling.  For  example,  the  oil  configuration  damping 
coefficient  provides  a  reasonable  estimate  for  future  gage  designs  that  use  an  FC-8  filled 
bellows.  Therefore  we  can  take  this  coefficient,  a  FEM  determined  vibration  frequency, 
and  make  predictions  of  new  gage  settling  times. 


Table  3-5.  CF2  vibration  frequency  calibration  results. 


(On  (rad/s) 

c 

Rise 

time(ms) 

Settling 
time  (ms) 

Baseline 

10996 

0.028 

0.17 

5.5 

Bellows 

11574 

0.027 

0.41 

2.2 

Oil 

7195 

0.043 

0.29 

8 

3.4.2  CF2  Virginia  Tech  Supersonic  Tunnel  Wall  Shear  Test 

The  CF2  gage  validation  testing  was  conducted  in  the  Virginia  Tech  Supersonic 
Tunnel.  Two  test  conditions  were  used:  a  pure  wall  shear  test  and  a  shock-impingement 
test.  The  total  pressure  set  point  was  1.03  MPa.  The  total  temperature  was  atmospheric. 
Total  pressure,  (P0),  wall  static  pressure  (Pst),  and  wall  static  temperature  Tst,  were 
measured  along  with  the  CF2  output. 

Each  CF2  Wheatstone  bridge  was  completed  with  a  Vishay®  Model  2310  signal 
conditioner.  These  amplifiers  were  set  using  the  10  Hz  low-pass  Butterworth  filter,  a  5  V 
input  voltage,  and  a  10k  gain.  A  National  Instruments  ®  data  acquisition  system  was 
controlled  by  programs  written  in  Lab  VIEW  ®.  The  sample  rate  was  500  Hz. 

Multiple  validation  tests  were  conducted  for  the  pure  wall  shear  case.  Figure  3-19 
shows  the  voltage  output  from  each  half-bridge  during  one  such  run.  The  red  traces 
indicate  output  from  the  head  half-bridge,  near  the  top  of  the  beam.  The  blue  lines 
indicate  traces  from  the  root  half-bridge,  near  the  bottom.  Dashed  lines  show  the 
expected  voltage  output  for  the  head  and  root  for  a  pure  85  Pa  wall  shear,  based  on 
calibration.  First,  let’s  follow  the  root  trace  with  time.  The  run  begins  without  flow, 
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indicated  by  zero  voltage  output.  After  5  s,  the  tunnel  starts  as  evidenced  by  the  jump  in 
voltage  by  the  root  bridge.  Near  17  s,  the  tunnel  pressure  starts  to  drop,  until  it  unstarts  at 
20s.  The  root  voltage  clearly  indicates  this  by  returning  to  zero.  Throughout  the  test,  the 
root  output  behaves  very  close  to  expectations. 


Vo  Root,  Vo  Head  vs.  Time 


Figure  3-19.  Pure  wall  shear  results  voltage  output  vs.  time. 

Turning  our  attention  to  the  red,  or  head  trace,  we  see  that  the  voltage  trends  do 
not  match  our  expectations.  In  fact,  this  output  trends  with  temperature.  This  is  best 
shown  in  Figure  3-20.  In  this  figure,  we  see  a  relationship  between  the  measured  static 
temperature  and  the  head  bridge  output  voltage.  Obviously,  the  temperature  is  a  source 
of  apparent  strain  for  these  test  conditions.  The  relatively  long  tunnel  start  and  test  times 
allow  for  heat  transfer  through  the  gage  head  and  down  the  beam.  The  first  order  time 
constant  comparing  the  peak  voltage  and  temperature  is  about  0.34  s.  With  the  long  test 
window,  the  gage  has  time  to  react  to  static  temperature  changes. 


44 


Temperature  and  Head  Bridge  Output  During  VT  SST  Mach=4  Run 


Figure  3-20.  Comparison  of  CF2  head  voltage  and  static  temperature  during  M4  run. 

Several  methods  of  post-process  temperature  correction  were  considered  and 
discarded.  Rather  than  develop  a  temperature  compensation  method,  which  would 
require  careful  temperature  measurements  and  physical  considerations  for  each  test  case, 
a  bounded  solution  method  is  developed.  This  method  considers  the  full  range  of 
possible  wall  shear  and  moment  solutions  comparing  test  and  calibration  data.  For  the 
test  case  presented  in  Figure  3-19,  we  are  only  interested  in  the  steady  state  conditions 
between  5  and  17  s.  This  window  of  interest  is  bound  by  the  plus  ‘+’  symbols.  Within 
this  region,  the  extrema  for  both  root  and  head  voltage  are  selected  as  bounds.  These 
extrema  are  indicated  by  four  circles.  The  solutions  to  Equations  2-3  and  2-4  provide 
corresponding  moment  and  wall  shear  values  for  each  of  these  circles.  These  solutions 
are  shown  in  Figure  3-21  as  points  A,  B,  C,  and  D. 


45 


Moment  (Nm) 

Figure  3-21.  CF2  pure  wall  shear  bound  solution. 


The  solutions  have  been  bound.  From  point  A  and  B,  the  wall  shear  lies  within 
the  range  of  6  to  220  Pa.  This  is  obviously  not  useful.  However,  since  we  know  the  head 
output  voltage  is  dominated  by  temperature  effects,  we  can  deduce  that  moment  is 
minimal.  With  this  knowledge,  we  can  move  to  the  point  of  the  curve  where  moment  is 
zero.  Examining  the  single  standard  deviation  solution  at  zero  moment  provides  a  very 
tight  estimation  for  wall  shear  indicated  by  point  K.  This  value  is  85  Pa  +  3  Pa.  This 
solution  fairly  well  matches  predicted  values  for  wall  shear  under  these  tunnel  conditions 
in  Table  1-3. 


3.4.3  CF2  Virginia  Tech  Supersonic  Tunnel  Shock  Impingement  Test 

A  wedge  validation  test  shadowgraph  is  superimposed  on  the  test  diagram  in 
Figure  3-22.  This  image  shows  how  an  oblique  shock  impinges  the  top  of  the  gage. 


Figure  3-22.  CF2  shock  validation  shadowgraph. 
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One  set  of  CF2  results  for  this  condition  is  shown  in  Figure  3-23.  The  root 
voltage  output  shown  in  blue  exhibits  similar  behavior  to  that  of  the  pure  wall  shear  case. 
Though  the  head  voltage  output  in  red  shows  some  trend  with  temperature,  the  signal  is 
now  dominated  by  moment. 


Figure  3-23.  CF2  shock  validation  voltage  output  versus  time. 


The  bounding  applied  to  these  output  signals  is  presented  in  Figure  3-24.  This 
time,  the  data  was  bound  to  the  one-sigma  solution.  For  example,  the  wall  shear  between 
points  E  and  F  ranged  from  -156  Pa  to  -13  Pa.  The  average  at  point  I  was  -84  Pa.  This 
negative  range  for  wall  shear  suggests  separation  as  might  be  expected  in  these 
conditions. 

A  significant  part  of  these  results  is  that  the  moment  solution  is  bound  between 
0.0002  Nm  and  0.0003  Nm.  Comparing  to  our  predicted  value,  this  measurement  is  of 
the  same  sign  and  magnitude.  For  the  first  time  using  this  new  gage  we  can  accurately 
state  that  the  gage  experienced  a  significant  pressure  differential  across  the  head.  Now 
combining  this  information  with  the  measured  negative  wall  shear,  the  experimentalist 
can  conclude  that  a  shock-impingement  has  caused  a  separation  of  the  flow.  No  previous 
skin  friction  gage  could  do  so.  In  fact,  earlier  direct  measure  gages  under  these  same 
flow  conditions  would  only  indicate  a  large  positive  wall  shear,  exactly  the  wrong 
conclusion  for  these  conditions.  Therefore,  this  new  gage  design  using  separate  half- 
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bridges  to  account  for  moment  has  proven  its  unique  utility  for  these  difficult  flow 
environments. 


Wedge  Mo  vs  Tau 
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Mo  (Nm) 
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A 

0.000105 

116.5052' 

B 

0.000354 

-299.994 

C 

0.000125 

-19.5779 

■D 

0.000334 

-163.911  ■ 

E 

0.000191 

-12.8782 

F 

0.000276 

-155.744 

G 

0.000198 

-61.2996 

H 

0.000268 

-107.322 

■1 

0.000233 

-84.3109- 

■'■Sir 
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Moment  i 


x  10' 


Figure  3-24.  CF2  bound  shock  impingement  solution. 
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4  Gage  Characterization 

After  CF2  testing,  we  recognized  the  need  to  develop  an  encompassing  design 
methodology.  Up  to  this  point,  the  designs  have  been  developed  through  piecemeal 
investigations.  This  chapter  builds  the  set  of  design  techniques  that  will  be  used  for 
optimization. 

4.1  Model  Results 

In  the  Introduction,  we  specified  a  set  of  nominal  load  conditions  in  Table  1-4. 
This  chapter  develops  the  expected  response  of  CF1  through  CF4  to  those  conditions  and 
establishes  the  fundamental  models  for  optimization.  Table  4-1  summarizes  the  results  of 
this  investigation.  In  the  table,  the  analysis  and  model  type  for  each  load  case  is 
presented.  For  each  gage  the  critical  output  parameters  are  listed.  We  will  now  examine 
how  each  of  these  parameters  was  calculated. 

Table  4-1.  CF1,  CF2,  CF3,  and  CF4  expected  output  to  nominal  loading. 


Analysis 

Model 

Load 

Output 

CF1 

CF2 

CF3 

CF4 

Static 

xw 

sx  (pa) 

13.31 

52.63 

11.85 

25.53 

3-D  FEM 

2000  Pa 

Sx  (gm) 

0.6 

6.1 

3.1 

9.5 

Static 

Xw 

St  (PS) 

49.3 

Calibrated 

2000  Pa 

Static 

M=10  NS 

M0  (mN-m) 

1.23 

4.16 

7.22 

7.22 

3-D  FEM 

Pl=500  Pa 

8Mo(ps) 

42.69 

143.17 

13.45 

26.34 

P3=58.25  kPa 

SMo(pm) 

2.6 

23.9 

6.0 

15.7 

Static 

M=10  NS 

8Mo(p£) 

132.7 

Calibrated 

Static 

Pw 

Sp  (pe) 

FAIL 

0 

0  Est. 

0  Est. 

Calibrated 

500  Pa 

Static 

Acc 

xa  (Pa) 

21.4 

20.6 

181.1 

75.0 

Solid 

lg 

Sa  (pc) 

0.1 

0.6 

0.8 

0.5 

Mechanics 

Xa/  Xw  (%) 

1.1 

1.0 

9.1 

3.7 

Thermal 

9V 

AT  (K) 

4 

3 

<1 

1 

3-D  FEM 

100  (kW/m2) 

sT  (pe) 

96 

72 

1.3 

24 

Dynamic 

f  (Hz) 

f  (Hz) 

3569.9 

1133.4 

1052.8 

1030.1 

Axi  FEM 

cd„  (rad/s) 

22430 

7121 

6615 

6472 

Est.  C 

0.027 

0.043 

0.043 

0.043 

Dynamic 

Step  Load 

Rise  (ms) 

0.2 

0.2 

0.2 

0.2 

Laplace 

Settle  (ms) 

3.6 

7.5 

7.9 

8.0 

Dynamic 

Step  Load 

Rise  (ms) 

0.29 

Calibrated 

Settle  (ms) 

8.0 
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The  response  of  each  gage  depends  on  its  geometry  and  material.  A  generalized 
sketch  of  the  fundamental  geometries  is  shown  in  Figure  4-1.  Table  4-2  lists  each  gage’s 
corresponding  fundamental  measurements.  Critical  structural  material  properties  are 
listed  in  Table  4-3.  Critical  thermal  material  properties  are  listed  in  Table  4-4.  A  quick 
reference  for  important  solid  mechanics  equations  is  provided  in  Table  4-5.  Finally,  the 
relevant  instrumentation  for  each  gage  is  shown  in  Table  4-6. 


Figure  4-1.  Geometric  parameters. 


Table  4-2.  Geometric  quick  reference. 


CF1 

CF2 

CF3 

CF4 

Material 

Aluminum 

Aluminum 

INVAR 

Aluminum 

Head 

rh 

3.175E-3  m 

4.7625E-3  m 

5.725E-3  m 

5.725E-3  m 

thh 

6.35E-4  m 

6.35E-4  m 

1 .27E-3  m 

1.27E-3  m 

Beam 

w 

2.54E-3  m 

2.54E-3  m 

3.048E-3  m 

3.048E-3  m 

th 

1.016E-3  m 

1.016E-3  m 

2.703E-3  m 

3.048E-3  m 

L 

6.858E-3  m 

1 .206E-2  m 

3.139E-2  m 

3.962E-2  m 

Table  4-3.  Structural  material  properties  reference.23 


Material  Properties 

Aluminum 

Nickel 

INVAR 

Tensile  Modulus  E  (Pa) 

7.00  E  10 

2.07  E  11 

1 .48  E  1 1 

Density  p  (kg/m3) 

2700 

8880 

8050 

Poisson's  Ratio  v  (m/m) 

0.33 

0.31 

0.23 

Coefficient  of  Thermal 
Expansion  CTE  (pm/mK) 

24 

N/A 

1.3 
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Table  4-4.  Thermal  material  properties  reference.24 


A1 

Nickel 

Invar 

Oil 

Air 

Thermal  Conductivity  k  (W/mK) 

210 

60.7 

10.15 

0.138 

0.026 

Specific  Heat  cp  (J/kgK) 

900 

460 

515 

1670 

1000 

Thermal  Diffusivity  a  (m2/s) 

8.6  E-5 

1.5  E-5 

2.4  E-6 

1.0  E-7 

2.23  E-5 

Table  4-5.  Solid  mechanics  equations.25 


Geometry 

Mass 

Equivalent  Acceleration  Shear 

Lf .  —  L  +  thh 

tot  n  «  , 

A  _„.r2  I  =  wth/\2 
~  71  rh  4 

Ac  =  w-th  3n  h 

mh  =  AhthhP 

mb  =AcLtotP 

Ta  =■ 

r  (r  th\ 

m  h  L tot  ^ 

V  1  J 

L 

+  mb  — 

b  2  J 

a 

L tot 

Table  4-6.  Gage  components  and  instrumentation.26’27 


CF1 

CF2 

CF3 

CF4 

Servometer  ®  Bellows 

FC-3-L 

FC-8-L 

FC-8-L 

FC-8-L 

Wheatstone  Bridge 

2  Half 

2  Half 

TBD 

TBD 

Vishay  Strain  Gage 

EA- 13-031  -CE-3  5  0 

E  A-  13-031-CE-350 

TBD 

TBD 

Rg 

350.0  ±  0.2% 

350.0  ±  0.2% 

TBD 

TBD 

GF 

2.11  ±  1.0% 

2.11  +  1.0% 

TBD 

TBD 

Vishay  2310  Gain 

10k 

10k 

TBD 

TBD 

Filter 

10  Hz  (VT) 

10  Hz  (VT) 

TBD 

TBD 

Vin 

5  V 

5  V 

TBD 

TBD 

4.2  Static  Response 

The  first  analysis  type  is  static.  In  this  case,  we  ignore  whether  the  load  is 
stepped,  ramped,  or  time-varying  and  only  treat  the  steady  state  case.  The  loads  of 
interest  are  wall  shear,  moment,  acceleration,  and  pressure.  Each  of  these  will  be 
addressed  in  turn.  The  analytic  cantilever  beam  solutions  provide  a  good  sanity  check  for 
more  complicated  models.  The  cantilever  beam  static  solutions  are  presented  in  Table 
4-7. 
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Table  4-7.  Cantilever  beam  static  load  solutions. 


Surface 

Load 

> 

> 

k 

ITT 

lg 

tw=2000  Pa 

Pi=500  Pa 

P3=  58250  Pa 

Pw= 

500  Pa 

Equation28 

Mr  =  *wAhLto, 

Mr  ■ th 
°T  -  21 

£  _  ^w^h^tot 

3EI 

M  =  2^~P^ 

°  3 

M0  ■ th 

<7mo  ~  21 

_  GMo 

Mo  E 

c  _  MoAot 

M°  2EJ 

C| 

II  II 

Wall  Shear 
Equations 

Output 

Si  (pS) 

ST  (pm) 

Smo (pe) 

8Mo  (pm) 

^Pw  (10  |as) 

Ta(Pa) 

Sa  (PS) 

CF1 

9.9 

0.4 

25.6 

1.4 

4.5 

22.9 

0.1 

CF2 

59.2 

6.3 

136.0 

21.6 

7.1 

21.9 

0.6 

CF3 

8.8 

2.3 

9.4 

3.7 

2.4 

193.6 

0.8 

CF4 

12.1 

4.4 

10.3 

5.7 

3.4 

79.0 

0.5 

4.2.1  Static  Wall  Shear 

The  first  load  for  consideration  is  the  nominal  2000  Pa  static  wall  shear.  Because 
we  are  interested  in  the  strain  distribution  on  the  surface,  a  full  3-D  finite  element  model 
(FEM)  is  required,  and  the  commercial  FEM  Code  ANSYS  ®  was  used  for  this  and  other 
purposes.29  We  will  use  CF2  as  a  case  study  since  calibration  data  is  available.  The 
finite  element  model  uses  the  Solid  92  element,  which  is  a  3-D,  10-node  structural 
tetrahedral  with  three  degrees  of  freedom,  i.e.  displacement  in  each  Cartesian  direction: 
ux,  uy,  and  uz.  The  beam  bottom  surface  nodes  were  held  to  zero  displacement.  The  wall 
shear  was  evenly  distributed  across  the  top  surface.  The  free  mesh  was  produced  using 
ANSYS  smart  meshing  set  to  6,  which  produced  1,366  nodes.  The  meshed  model  and 
nodal  solution  are  shown  in  Figure  4-2. 
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Figure  4-2.  CF2  3-D  FEM  result  for  2000  Pa  wall  shear. 


Table  4-8  presents  the  Voltage-Pascal  ratios  for  the  root  and  head  from  static  wall 
shear  calibration.  These  are  the  quantities  cn  and  C21  for  Equations  2-3  and  2-4. 
Assuming  ideal  instrumentation,  we  can  back  out  the  expected  microstrain  for  this  load 
using  Equation  3- 1 .  The  manufacturer  published  gage  factor  is  2. 1 1 ,  and  the  23 1 0  Signal 
conditioner  is  set  with  Gain=10,000  and  Vin=5  V.  From  this  we  expect  to  see  a  strain  of 
49.3  ps  at  the  root,  positive  for  the  tension  side  and  negative  for  compression  side. 

2-V0 

S  = -  Equation  4-1 

Vin  GF  ■  Gain 


We  expected  the  maximum  strain  calculated  using  ANSYS  to  be  higher  than  from 
the  actual  gage.  This  is  because  each  strain  gage  measures  the  strain  over  a  relative  large 
area  compared  to  the  point  of  maximum  strain.  The  fact  that  our  model  is  within  10%  of 
the  measured  strain  is  good.  It  shows  that  strain  gages  were  well  placed  during 
fabrication.  It  also  gives  confidence  in  the  modeling  technique. 


Table  4-8.  Expected  CF2  output  to  200  Pa  wall  shear. 


Facility 

Root 

Head 

Variation  from  Ideal 

Calibration  Factors 

0.0013  V/Pa 

0.0003  V/Pa 

Output  for  2000  Pa 

2.6  V 

0.6  V 

Ideal  Bridge 

49.3  ps 

1 1 .4  ps 

0% 

Cantilever  Beam 

59.2  ps 

N/A 

20% 

3-D  FEM 

52.6  ps 

N/A 

6.7  % 
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4.2.2  Static  Moment 


The  next  load  for  consideration  is  a  severe  Mach  10  normal  shock  impinging  on 
the  center  of  the  gage  head  which  for  CF2  yields  a  0.00416  N-m  moment.  Again, 
because  we  are  interested  in  the  strain  distribution  on  the  strain  gage  attachment  surfaces, 
we  have  to  use  a  full  3-D  finite  element  model.  We  continue  with  the  CF2  as  a  case 
study,  since  calibration  data  is  available.  The  finite  element  model  uses  the  same  Solid 
92  element.  The  beam  bottom  surface  nodes  were  held  to  zero  displacement.  The 
moment  load  was  applied  as  an  axially  varying  pressure  load  distributed  across  the  top 
surface.  The  free  mesh  was  produced  using  ANSYS  smart  meshing  set  to  6,  which 
produced  1,366  nodes.  The  meshed  model  and  nodal  solution  are  shown  in  Figure  4-3. 

-219.0  — | 

-179.0  — 

-138.0  — 

-98.2  — 

-58.0  — 

-17.7  — 

22.5  — 

62.7  _ 

103.0  _ 

143.0  ^ 

Figure  4-3.  CF2  3-D  FEM  result  for  impinging  Mach  10  normal  shock. 


Table  4-9.  Expected  CF2  output  to  a  Mach  10  normal  shock. 


Facility 

Root 

Head 

Variation 

Calibration  Factors 

1677  V/N-m 

1806  V/N-m 

Output  due  to  Ml 0  NS 

7.0  V 

7.5  V 

Ideal  Bridge 

132.7  ps 

142.2  ps 

0% 

Cantilever  Beam 

136.0  ps 

2.5  % 

3-D  FEM 

143.17  ps 

7.9  % 

Table  4-9  presents  the  resulting  Voltage-Pascal  ratios  for  the  root  and  head  from 
static  moment  calibration.  These  quantities  are  C12  and  C22  from  Equations  2-3  and  2-4. 
Assuming  ideal  instrumentation,  we  can  back  out  the  expected  microstrain  for  this  load 
using  Equation  3-1  with  previous  settings.  From  this  we  expect  to  see  a  strain  of  132.7 
ps  at  the  root,  positive  for  the  tension  side  and  negative  for  compression  side.  The 
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ANSYS  calculated  maximum  strain  should  be  higher  than  from  the  actual  gage.  Again 
the  strain  gage  covers  more  area.  Even  so,  the  model  is  still  within  10%  of  the  measured 
strain. 

4.2.3  Inertia  Loading 

The  last  load  for  consideration  is  inertia.  This  is  accomplished  by  calculating  an 
equivalent  acceleration  shear  due  to  inertia,  xa.  The  development  of  this  calculation  is 
shown  graphically  in  Figure  4-4.  In  this  equation,  a  is  the  measure  of  acceleration  in 
m/s2.  Our  nominal  maximum  acceleration  is  1-g  or  9.8  m/s2.  By  applying  xa  to  the  beam 
surface,  we  achieve  the  same  strain  output  near  the  base  of  the  beam  as  caused  by  inertia. 
Because  the  gage  is  performing  within  the  elastic  region  we  can  simply  apply  the  ratio 
Ta/xw  to  the  strain  output  for  the  pure  wall  shear  case.  As  an  example,  should  CF2 
experience  a  1-g  acceleration  and  2000  Pa  wall  shear,  it  would  see  a  relative  strain  due  to 
inertia  of  2 1 .4/2000  ex,  or  about  0.1  ps. 


Figure  4-4.  Calculation  for  equivalent  inertial  shear. 
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4.2.4  Static  Pressure 

The  gage  design  using  incompressible  fluid  filler  in  the  bellows  eliminates 
uniform  wall  pressure  as  a  source  of  strain.  A  vacuum  test  calibration  of  the  oil-filled 
CF2  proved  that  there  is  no  voltage  output  due  to  changing  uniform  pressure.  Therefore, 
for  our  model  considerations,  uniform  pressure  loading  on  CF2,  CF3,  and  CF4  is 
assumed  to  produce  zero  microstrain  output.  The  CF1  case  failed  to  eliminate  uniform 
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wall  pressure  as  a  variable.  Therefore,  pressure  is  an  indeterminate  source  of  apparent 
strain  for  CF 1 . 

4.3  Dynamic  Response 

When  you  first  step  on  a  scale  in  the  bathroom,  it  bounces  around  and  finally  settles 
on  your  weight.  This  under-damped  second  order  system  response  is  similar  to  what 
happens  with  these  direct  measure  gages.  They  are  under-damped,  which  means  they 
oscillate  with  an  exponential  decay  toward  a  steady  load  value.  In  order  to  characterize 
how  these  gages  respond,  we  really  want  to  know  how  long  the  gage  takes  to  get  within 
+  10%  of  the  steady  value.  This  is  called  the  gage  settling  time.  In  order  to  determine 
this  key  gage  parameter,  we  must  look  at  two  aspects  of  dynamic  response.  First,  we 
need  to  characterize  each  gage’s  dynamic  coefficients:  natural  frequency,  oon,  and 
damping  ratio,  C,.  Second  we  have  to  develop  solutions  to  the  ordinary  differential 
equations  for  our  different  expected  loads  cast  as  forcing  functions.30 

4.3.1  Vibrational  Frequency 

The  first  model  for  verification  is  1st  mode  bending  response.  We  continue  the 
trend  of  employing  CF2  due  to  its  available  calibration  data.  However,  in  this  case,  we 
use  a  poor  approximation  of  just  the  gage  beam  without  the  head.  This  way  we  can  find 
some  useful  relations  for  finite  element  modeling.  The  natural  frequency  of  just  the  beam 
is  calculated  using  the  analytic  solution  for  a  clamp-free  simple  cantilever  beam.  Table 
4-10  shows  the  clamp-free  response  of  this  simple  cantilever  beam.  Additionally,  the 
table  provides  a  means  to  calculate  dimensions  for  a  cylindrical  beam  of  equivalent 
natural  frequency.  Now  we  can  test  whether  this  equivalent  cylindrical  beam  provides 
reasonable  output  for  a  variety  of  finite  element  models. 

Several  finite  element  models  are  developed  to  match  the  proposed  rectangular 
aluminum  beam.  These  are  presented  in  Table  4-11.  Clearly,  these  models  well  match 
the  theoretical  result.  All  of  the  models  are  within  2%  of  the  theoretical  result  which 
verifies  the  coding  technique  applied  in  ANSYS.  The  axisymmetric  case  checked  well 
against  the  3-D  circular  model.  The  axisymmetric  model  is  important  because  it  provides 
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an  accurate  solution  for  vibration  frequency  with  a  greatly  reduced  number  of  elements 
compared  to  the  full  3-D  case. 


Table  4-10.  Cantilever  beam  solution  for  1st  mode  bending  response. 


Governing  Equation 
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12  -w-th 


n  ■  r 


4 n  ■  r2 


eg 


=  -Jl/3th 


Output 
f=  5182.1  Hz 
can  =  32,560  rad/s 
req=5.87E-4m 


Table  4-11.  Comparison  of  several  FEMs  to  cantilever  theory  for  1st  mode  bending. 


1-D  FEM 


Axi-Harmonic 


3-D  Cylinder 


3-D  Block 


ANSYS 

Element 

Type 


Beam 

Theory 


Beam3 


Plane83 


Solidl86 
Mopt  Pyra 


Solidl86 


Model 


X 


T 

L 


-tot 


f(Hz) 


5,182.1 


5,175.3 


5183.9 


5187.5 


5248.6 


on  (rad/s) 


32,560 


32,517 


32,571 


32,594 


32,978 


Variation 


0.1% 


-0.0% 


-0.1% 


-1.3% 


Now  that  we  know  the  axisymmetric  finite  element  model  is  sufficient  to 
determine  natural  frequency,  we  can  develop  accurate,  but  computationally  small  models 
for  the  real  gages.  The  axisymmetric  model  allows  us  to  model  the  full  beam  and 
bellows.  The  ANSYS  element  used  for  this  analysis  was  the  Plane  83.  This  axiharmonic 
element  is  simply  an  axisymmetric  element  on  which  non-axisymmetric  loads  may  be 
applied.  For  example,  to  model  a  surface  wall  shear,  the  mode  is  set  to  1,1  and  the  load  is 
applied  as  fx=  xw  and  fz=-Tw.  The  meshed  model  and  first  mode  bending  shape  are  shown 
in 
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Table  4-12.  The  measured  vibration  frequency  for  the  CF2  gage  is  included.  We  can  see 
that  the  finite  element  model  performs  well,  predicting  the  response  to  1.0  %.  The 
damping  coefficient  for  CF2  comes  directly  from  measurement  Figure  3-18.  This  data 
was  used  to  make  assumptions  about  the  damping  coefficient  for  each  of  the  other  gages. 
CF1  has  lower  damping  since  it  lacks  oil- fill. 


Table  4-12.  Natural  frequency  output  for  each  gage. 


CF1 

CF2 

CF3 

CF4 

Axiharmonic 

FEM 

& 

,\ 

In 

& 

5 

- 

i 

if 

w 

I 

-1 

Node  count 

1215 

1982 

1922 

1925 

c 

0.027 

0.043 

0.043 

0.043 

f(Hz) 

3569.9 

1133.4 

1052.8 

1030.1 

con  (rad/s) 

22,430 

7121 

6615 

6472 

Calibration 

7195  rad/s 

Variation 

1.0% 

4.3.2  Step  Load 

We’ve  already  stated  that  this  gage  exhibits  second  order  system  response. 
Equation  4-2  is  the  governing  2nd  order  linear  non-homogeneous  ODE.  This  equation 
shows  how  the  system  will  move  in  the  x-direction  given  a  time-varying  forcing  function, 
F(t),  and  two  separate  initial  conditions.  If  we  divide  through  by  mass,  we  get  a  familiar 
form  with  c=2^(on  and  k=ton  ,  where  C,  is  the  damping  factor  and  oon  is  the  circular 

33 

frequency  in  rad/s.  In  order  to  verify  our  models,  we  use  con  =  7195  rad/s  and  C,  =  0.043 
for  the  CF2  gage.  For  a  step  wall  shear  input,  a  Laplace  transform  time  domain  solution  is 
shown  in  Equation  4-4. 

Mx  +  Cx  +  Kx  =  F(t )  Equation  4-234 

x  +  cx  +  kx  =  f{t )  Equation  4-3 
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x(t )  =  u  —  ue  ^ 


f 


Equation  4-435 


The  transform  solution  is  compared  with  the  calibration  data  for  the  oil-filled  CF2 
in  Figure  4-5.  The  transformation  solution  overshoots  the  peak  value  achieved  during 
experiment.  Since  we  are  using  a  transformation  model  for  optimization,  this  allows  a 
conservative  estimate  for  the  maximum  displacement.  In  other  words,  we  will 
automatically  design  a  gap  size  that  will  accommodate  the  overshoot  without  railing. 


Oil 


Figure  4-5.  CF2  step  load  response.  Comparison  between  calibration  data  and  transformation 
solutions. 


Unfortunately,  the  model  under  predicts  the  signal  rise  time,  but  more  critically, 
the  settling  time.  This  must  be  accounted  for  when  choosing  a  minimum  vibration 
frequency  for  settling  time.  The  full  parameters  are  presented  in  Table  4-13. 
Conservatively,  we  can  apply  an  85%  design  factor  to  the  total  run  time  to  ensure  the 
gage  settles  as  desired. 

Now  that  we  have  a  good  method  to  determine  the  rise  and  settling  time 
using  the  Laplace  transform  solution,  we  can  characterize  the  dynamic  response  of  each 
gage.  Figure  4-6  shows  the  expected  dynamic  response  for  each  of  our  gages  using  the 
coefficients  from  Table  4-12.  Clearly  the  CF1  gage  settles  the  quickest  just  from  quick 
observation.  The  numerical  results  for  rise  and  settling  time  for  each  gage  are  listed  as 
well.  Clearly  all  gages  will  settle  sufficiently  to  make  an  accurate  measurement  within 
10  ms. 
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Table  4-13.  CF2  Validation  data. 


Experiment 

Laplace 

Variation 

Safety  Factor 

Rise  Time 

0.000299 

0.000208 

30.4  % 

Settling  Time 

0.008268 

0.007059 

14.6  % 

0.85  (min  test  time) 

Peak  Value 

1.535 

1.872 

-22  % 

2  (max  tw  static) 

Figure  4-6.  Dynamic  response  to  step  load. 


4.4  Thermal  Response 

The  thermal  response  requires  a  transient  analysis.  Because  of  the  interaction  of 
the  gage,  oil,  and  air,  a  3-D  finite  element  model  is  used.  Additionally,  because  of  the 
effect  of  surrounding  material,  such  as  a  floor  plate,  a  larger  model  is  developed  for  each 
gage.  We  assume  the  heat  flux  load  of  100  kW/m  acts  on  the  gage  prior  to  achieving 
steady  state  during  a  test  run.  We  must  account  for  these  transient  effects,  such  as  tunnel 
start  for  example.  This  is  done  by  applying  the  heat  flux  for  a  full  0. 1  s  instead  of  just  the 
10  ms  steady  state  time. 
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The  ANSYS  element  type  for  this  model  is  the  Solid  98,  which  is  a  3-D,  10-node 
tetrahedral.  The  element  is  set  to  only  allow  one  degree  of  freedom,  temperature.  As  an 
additional  check  on  potential  thermal  gradients,  the  heat  flux  is  varied  axially  across  the 
surface  of  the  model.  The  nominal  load  is  applied  at  axial  centerline.  The  leading  edge 
is  2.5%  greater,  and  the  trailing  edge  is  2.5%  lower.  This  variation  of  5%  across  the 
model  surface  mimics  boundary  layer  growth  and  allows  us  to  investigate  temperature 
variations  through  the  beam  thickness.  Figure  4-7  shows  the  applied  spatially  varying 
heat  flux  loads  for  each  model. 


□ 

□ 

□□ 

□□ 

□ 

□□ 


Figure  4-7.  5%  spatially  varying  heat  flux,  nominal  100  kW/m2. 


The  top  left  image  shows  the  CF1  model.  The  smart  mesh  option  of  10  yielded 
3512  elements.  The  top  right  image  shows  the  CF2  model.  The  smart  mesh  10  option 
yielded  3906  elements.  The  bottom  left  model  shows  the  CF3  gage  in  an  INVAR  plate. 
The  smart  mesh  9  option  for  this  model  produced  9368  elements.  Finally,  the  bottom 
right  image  shows  the  CF4  gage.  The  smart  mesh  9  option  produced  14,734  elements. 

All  models  had  an  initial  temperature  of  289  K. 

Performing  the  transient  analysis  for  0.1  s  produced  the  thermal  distributions 
shown  in  Figure  4-8.  The  arrows  for  CF1  and  CF2  point  to  the  approximate  vertical 
location  of  the  head  strain  gages.  We  can  see  that  this  particular  set  of  gages  experience  a 
relatively  high  amount  of  temperature  gradient  due  to  this  nominal  load.  Looking  at  the 
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CF3  and  CF4  gages,  we  see  large  regions  where  little  temperature  change  occurs.  Strain 
gages  placed  within  this  region  should  experience  little  apparent  thermal  strain. 


Figure  4-8.  Temperature  distribution  100  kW/m2  for  0.1  s. 


In  order  to  quantify  the  amount  of  strain  experienced  at  these  locations,  the  simple 
formula  sT  =  CTE  •  AT  is  applied.36  The  small  thermal  expansion  coefficient  for 
INVAR  makes  a  large  impact  on  expected  thermal  strain  due  to  temperature.  The 
compiled  results  for  the  expected  response  at  the  ‘head’  gage  are  supplied  in  Table  4-14. 
The  term  ‘head’  is  placed  in  quotes  to  emphasize  the  new  lower  location  on  CF3  and 
CF4. 

Please  recall  that  if  both  sides  of  the  beam  see  the  same  strain  output,  we  should 
get  a  net  voltage  output  of  0V.  However,  for  arguments  sake,  if  the  beam  temperature 
varies  5%  from  front  to  back,  this  produces  a  measured  5%  As  across  the  bridge,  and  it 
produces  a  proportional  voltage  output.  In  order  to  grasp  that  effect,  the  amount  of 
thermal  strain  produced  from  a  5%  variation  is  also  shown  in  Table  4-14. 


Table  4-14.  Expected  strain  output  at  ‘head’  due  to  nominal  heat  flux. 


CF1 

CF2 

CF3 

CF4 

AT 

4K 

3K 

<1K 

<1K 

St 

96  ps 

72  ps 

1.3  ps 

24  ps 

Assumed  5%  Ast 

4.8  ps 

3.6  ps 

0.1  ps 

1.2  ps 
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4.5  Combined  Response 

The  measure  of  success  for  CF1  through  CF4  is  how  well  each  meets 
specifications.  The  design  load  conditions  are  presented  here  again  in  Table  4-15.  This 
time  we’ve  added  the  additional  complexity  of  spatially  varying  heat  flux  examined  in 
the  last  section.  Given  these  load  conditions,  how  close  does  each  of  these  gages  come  to 
measuring  the  2000  Pa  wall  shear? 

Table  4-15.  Nominal  design  loads 


Time 

(ms) 

T 

wx 

(Pa) 

°x 

(g) 

VvW 

(kW/m2) 

Load 

10 

2000 

1 

100 

Figure  4-9  shows  the  expected  magnitude  of  strain  from  these  nominal  loads 
assuming  a  5%  thermal  strain  output.  In  general  we  see  that  in  all  cases,  the  signal  is 
dominated  by  shear  strain,  which  is  a  good  thing.  Though  CF1  has  been  removed  from 
contention  due  to  pressure  effects,  we  see  that  it  would  have  suffered  due  to  thermal 
strain.  The  other  gages  have  relatively  small  amounts  of  thermal  strain.  As  expected  the 
INVAR,  CF3  gage,  performs  exceptionally  in  limiting  apparent  strain  due  to  thermal 
gradient.  CF2-CF4  experience  relatively  small  amounts  of  apparent  strain.  Even  in  the 
case  of  CF3,  which  has  a  proportionately  large  amount  of  inertial  strain,  it  produces  less 
than  10%  of  the  total  strain  output. 


Strain  Output 


□  Shear  ■  Inertia  □  Thermal  Gradient 


Figure  4-9.  Strain  magnitudes  due  to  nominal  wall  shear. 
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In  order  to  quantify  these  sources  of  strain,  we  look  at  the  relative  magnitude 
percentage  of  strain  over  the  strain  from  just  wall  shear.  This  is  shown  in  Table  4-16. 
This  table  shows  that  CF2,  CF3,  and  CF4  should  have  reasonably  small  apparent  strain 
under  our  nominal  load.  CF1  clearly  would  have  had  significant  thermal  issues.  For  the 
optimized  gages,  apparent  strain  provides  less  than  10  %  of  signal  output. 


Table  4-16.  Relative  magnitude  of  strain  output. 


l^a /  £t| 

|st/8x| 

CO 

Sr 

o 

£i/stot 

CF1 

i% 

36% 

i% 

26% 

CF2 

i% 

7% 

i% 

6% 

CF3 

9% 

1% 

8% 

1% 

CF4 

4% 

5% 

3% 

4% 

In  order  to  understand  how  these  sources  of  strain  affect  the  output  signal,  we 
have  to  consider  whether  the  strain  is  positive  or  negative.  In  the  case  of  inertia,  we 
know  that  it  opposes  the  motion  induced  by  wall  shear.  Therefore  for  simplicity  if  we 
consider  the  tension  side  of  the  beam,  we  expect  to  have  a  positive  strain  from  wall  shear 
and  a  negative  strain  from  inertia.  The  apparent  thermal  strain  could  be  positive  or 
negative.  Now  we  can  band  our  expected  output  signal  for  each  gage  due  to  the  nominal 
load  conditions.  This  is  done  in  Figure  4-10. 

Expected  Gage  Output 


CF1  CF2  CF3  CF4 

Figure  4-10.  Expected  output  for  nominal  conditions. 
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The  lighter  gages,  CF1  and  CF2,  have  an  uncertainty  band  that  is  spread  almost 
equally  about  the  true  value.  The  heavier  gages  are  likely  to  indicate  a  lower  value  of 
wall  shear  than  is  actually  measured  due  to  inertia.  The  expected  output  values  are  in 
Table  4-17. 


Table  4-17.  Output  measurements  to  nominal  load  case. 


Low  (Pa) 

2000  -  Low 

High  (Pa) 

1 

o 

o 

o 

<N 

Mean  (Pa) 

2000 

2000 

CF1 

1257 

37% 

2700 

-35  % 

1979 

CF2 

1843 

8% 

2116 

-6% 

1979 

CF3 

1808 

10% 

1830 

9% 

1819 

CF4 

1831 

8% 

2019 

-1  % 

1925 

An  experimentalist  using  one  of  these  gages  can  limit  uncertainty  due  to 
acceleration.  This  is  accomplished  by  measuring  acceleration  during  testing  of  the  base 
model  with  an  accelerometer.  The  most  conservative  approach  is  to  multiply  the 
maximum  measured  acceleration  by  the  coefficient  in  Table  4-18.  This  provides  the 
output  (Pa)  that  should  be  added  to  the  wall  shear  measured.  As  an  example  using  CF3,  a 
maximum  1  g  acceleration  of  9.8  m/s  is  multiplied  by  the  19.74  kg/m  acceleration 
coefficient.  This  yields  a  193.6  Pa  correction  that  should  be  added  to  the  raw  wall  shear 
measurement. 


Table  4-18.  Accelerometer  employment  with  direct  measure  gages. 


CF1 

CF2 

CF3 

CF4 

Acceleration  Coefficient  (kg/m  ) 

2.33 

2.23 

19.74 

8.05 

lg  equivalent  shear  (Pa) 

22.9 

21.9 

193.6 

79.0 

4.6  Gage  Characterization  Techniques 

This  section  summarizes  the  selected  gage  characterization  techniques.  These  are 
the  simplest  available  modeling  techniques  that  accurately  predict  gage  response.  In  the 
case  of  displacement,  simpler  techniques  could  suffice,  except  this  parameter  is  already 
available  after  examining  the  full  3-D  FEM  for  corresponding  strain.  Table  4-19  shows 
the  selected  model  techniques,  output,  and  calibration  method  for  verification. 
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This  table  captures  the  critical  parameters  for  direct  skin  friction  gage  design. 

The  five  sources  of  strain  are  each  examined:  xw,  M0,  AT,  a,  and  Pw.  Wall  pressure  is  not 
modeled  since  through  vacuum  calibration,  it  has  been  eliminated  as  a  source  of  output 
voltage.  The  final  parameter  to  be  determined  for  gage  performance  is  settling  time.  If 
we  assume  a  damping  coefficient,  we  can  estimate  the  gage  dynamic  response  by 
applying  a  FEM  predicted  vibrational  frequency. 

The  output  parameters  of  interest  are  listed.  There  are  five  strains  of  interest:  ex, 
£Mo,  sa,  Sat,  and  sPw.  As  with  wall  pressure,  the  corresponding  strain  is  neglected. 
Normally,  we  would  look  at  the  corresponding  displacements  for  each  of  these  strains. 
However,  we  are  only  interested  in  displacement  in  terms  of  impact  on  gap  size. 
Therefore,  we  ignore  pressure  and  thermal  induced  displacements  since  they  do  not  act  in 
a  way  to  impact  gap  sizing.  Acceleration  displacement  is  not  considered  for  two  reasons. 
First,  wall  shear  effect  on  gap  size  is  much  greater  than  acceleration  by  design.  This  is 
because  we  want  small  relative  displacement  due  to  acceleration.  Second,  acceleration 
opposes  wall  shear,  so  it  actually  decreases  the  necessary  gap.  Thus  our  displacements  of 
interest  have  been  narrowed  to  just  two:  Sx  and  8m0- 


Table  4-19.  Gage  characterization  summary. 


Input 

Tw 

Mo 

1st  mode 
excitation 

Ta 

Pw 

Model 

T 

3-D  FEM 

3-D  FEM 

3-D  FEM 

N/A 

Axi  FEM 
Estimate 

A 

i 

lec 

hanic 

;s 

Output 

£t 

8X 

SMo 

Smo 

con 

C  (Est.) 

Sa 

Sat 

0  V 

Verifi¬ 

cation 

¥ 

x|. 

Trt 

Dfrj 

The  last  output  parameters  considered  are  necessary  to  determine  the  gage  settling 
time.  Since  we  are  designing  within  a  family  of  gages,  all  using  oil-filled  bellows  and  of 
similar  dimensions,  an  estimate  for  damping  coefficient  is  made  from  CF2  dynamic 
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calibration  data.  This  estimated  constant  of  <^=0.043  actually  requires  no  further 
modeling.  With  this  parameter  established,  we  only  need  to  determine  expected 
vibrational  frequency,  con. 

Finally,  the  table  shows  the  verification  techniques  for  each  model.  The  static 
wall  shear,  moment,  and  natural  frequency  calibrations  showed  good  agreement  with 
FEM  for  CF2.  The  vacuum  test  showed  that  wall  pressure  changes  did  not  produce 
voltage  output  for  CF2.  Static  temperature  measurements  during  testing  ensure  correct 
heat  flux  modeling.  This  point  will  be  thoroughly  investigated  in  Chapter  6:  Simulations 
and  also  shows  good  agreement  for  CF2  testing.  The  equivalent  acceleration  can  be 
verified  through  a  combination  of  static  wall  shear  calibration  and  accelerometer 
measuring.  For  design  purposes,  we  assume  the  solid  mechanics  model  is  accurate. 
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5  Optimization 

During  each  design,  the  gage  geometry  and  material  properties  were 
parametrically  varied  within  constraints  to  achieve  desired  performance.  The  geometry 
and  material  properties  made  up  the  design  variables.  The  output  variables  such  as 
strains,  frequencies,  and  displacements  were  state  variables  that  could  be  cast  as  objective 
functions.  An  objective  function  is  designed  to  be  minimized.  As  an  example,  we  would 
like  to  minimize  the  difference  between  our  maximum  strain  output  and  the  proportional 
strain  limit  while  staying  below  this  limit.  Over  several  gage  iterations,  from  CF1 
through  CF4,  the  number  of  state  variables  considered  increased  as  well  as  the 
sophistication  of  the  design  approach.  This  reflected  the  understanding  of  the  physical 
complexity  of  the  environmental  interaction  of  the  flow,  gage,  and  measurement  devices. 
It  also  illustrated  that  load  conditions  substantially  change  the  design  of  direct  measure 
gages.  The  optimization  applied  in  developing  CF3  and  CF4  fully  constrained  the  design 
space  and  represents  a  true  optimization  scheme  considering  all  design  and  state 
variables. 

The  robustness  of  the  optimization  approach  developed  in  this  chapter  allows  it  to 
be  applied  to  several  design  scenarios.  Some  nominal  conditions  for  three  very  different 
test  conditions  are  presented  in  Table  5-1.  Here  we  see  conditions  for  a  high-enthalpy 
short  duration  test  condition  such  as  LENS.  Next  come  the  longer,  lower  wall  shear 
blow-down  facility  conditions.  Finally,  the  design  conditions  for  a  typical  flight  test  are 
presented.  The  optimization  approach  will  be  applied  to  each  of  these  conditions. 


Table  5-1  Test  regimes  for  gage  optimization. 


Test  Conditions 

Settle  time 

F 

(Hz) 

Tw 

(Pa) 

ax 

m/s2 

kW/m2 

Mo 

Nm 

Shock  (High  Enthalpy) 

10  ms 

1000 

2000 

9.8 

100 

0.00722 

Blow  Down 

1  s 

10 

85 

0.098 

-5 

0.00118 

Flight 

1  s 

10 

2000 

9.8 

100 

0.00118 
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5. 1  Optimization  for  Short  Duration  High  Enthalpy  Testing 

The  first  test  condition  to  consider  is  the  familiar  nominal  condition  investigated 
in  Chapter  4:  Gage  Characterization.  One  added  complication  is  the  possibility  of  trying 
to  measure  a  lower  wall  shear  of  1000  Pa.  The  optimization  approach  will  be  applied  to 
develop  an  INVAR  design,  CF3,  and  an  aluminum  design,  CF4. 

5.1.1  CF3  Gage  Design 

After  achieving  a  gage  that  detects  shock  impingement,  CF2,  there  was  a  desire  to 
greatly  reduce  the  impact  of  thermal  output  from  the  upper  half-bridge  in  the  pure  wall 
shear  testing  in  the  VT  supersonic  tunnel.  As  discussed  in  the  modeling  section,  the 
thermal-strain  effect  was  considered  a  possible  source  of  error.  By  turning  our  attention 
to  a  material  with  low  thermal  expansion  characteristics,  we  hoped  to  greatly  reduce  or 
eliminate  this  effect.  Figure  5-1  shows  some  various  configurations  and  the  resulting 
temperature  distribution  for  expected  LENS  heat  flux.  Although  polyethersulfone  was 
considered,  the  requirement  to  switch  to  semi-conductor  strain  gages  made  this  material  a 
poor  candidate.  Thermal  strain  is  hard  enough  to  deal  with  in  the  linear  regime.  Semi¬ 
conductor  strain  gages  have  expansion  rates  that  are  non-linearly  related  to  the  substrate. 
Therefore,  we  chose  INVAR  which  has  a  low  coefficient  of  thermal  expansion. 


Baseline 

MacorCap 

Cu  Heat  Sink 

Invar 

Aluminum 

0.0003m 

1cm 

'300  '  '302  '  '304  '306  700 

301  303  305  650 


Figure  5-1.  Thermal  analysis  for  temperature  gradient  at  LENS.  FEM  using  100  kW/m2  for  0.1  s. 

Along  with  the  material  change,  an  oil  leak  problem  required  a  geometric  change. 
Initially,  we  considered  TIG  welding  or  soldering  the  bellows  to  the  head  and  base  of  the 
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gage.  After  further  consultation  with  our  machinists,  we  decided  to  enlarge  the  head  and 
add  a  groove  to  seat  the  bellows  in  epoxy  at  the  head  and  base.  This  increased  the  gage 
head  thickness  and  diameter.  This  obviously  penalizes  us  in  terms  of  allowable  inertial 
loading,  but  solving  the  thermal-mechanical  problem  has  taken  priority  in  this  design. 

Finally,  we  wanted  the  flexibility  to  use  a  full  Wheatstone  bridge  configuration  at 
both  the  “head”  and  root  locations.  As  shown  in  the  theoretical  development,  this  allows 
for  further  temperature  compensation  at  the  “head”  and  root  of  the  beam.  Even  using  the 
smallest  available  commercial  gages  Vishay  ®  EA-00-015LA-120,  we  still  limited  the 
minimum  design  width  to  0. 12  in.  Notice  that  the  “head”  location  has  been  moved 
considerably  further  down  the  beam  away  from  the  flow  in  Figure  2-8.  This  is  done  to 
avoid  thermal  gradients,  hence  the  quotes. 

So  in  general  the  new  design  has  many  features  to  address  the  thermal-strain  and 
leaking  problems,  but  has  greatly  increased  the  overall  size  and  weight  of  the  gage  and 
raised  the  CG.  Table  5-2  shows  the  combined  impact  of  the  design  changes  prior  to 
optimization.  First,  we  see  that  the  increased  density  will  impact  the  head  and  beam 
weight.  The  beam  width  puts  a  20%  penalty  on  weight.  The  head  of  CF3  is  nearly  9 
times  heavier  than  that  of  CF2.  This  also  pushes  the  center  of  gravity  higher.  All  of 
these  design  changes  have  a  detrimental  impact  on  acceptable  inertial  loading.  All  of 
these  variations  show  the  impact  of  designing  with  thermal-strain  considerations  at  the 
forefront. 


Table  5-2.  Design  feature  impact  on  size.  Comparison  between  CF2  and  CF3. 


CF2 

CF3 

Change 

Density  (kg/m  ) 

2700 

8050 

3  x  increase 

Beam  width  (in) 

0.1 

0.12 

1 ,2x  increase 

Head  thickness  (in) 

0.025 

0.05 

2  x  increase 

Head  diameter  (in) 

0.375 

0.4508 

1 ,2x  increase 

Head  volume  (m3) 

4.53E-8 

1.31E-7 

2.9x  increase 

Head  mass  (kg) 

1.22E-4 

1.05E-3 

8.6x  increase 

Despite  the  acceleration  limitation,  we  can  still  design  a  CF3  gage  that  achieves 
all  of  our  test  design  objectives.  In  order  to  determine  those  objectives,  we  need  to  settle 
on  reasonable  test  conditions.  This  gage  design  is  for  LENS  facilities.  Unfortunately,  the 
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variation  in  flow  conditions  across  these  facilities  is  rather  extreme.  These  conditions 
considered  are  presented  in  Table  5-3.  After  discussion  with  the  LENS  personnel,  we 
decided  to  build  a  gage  that  could  meet  all  of  the  LENS  I  and  LENS  II  requirements. 

This  meant  being  able  to  settle  within  10  ms  and  measure  a  1000  Pa  wall  shear,  without 
having  problems  under  the  higher  load  condition.  This  will  be  examined  in  considerable 
detail  as  we  examine  each  state  variable.  This  effectively  means  we  are  designing  for 
three  loading  conditions:  Case  1)  xwmin=1000  Pa  Case  2)  Twmax=2000  Pa  and  Case  3) 
M=10  shock  impingement.  Although  extreme,  the  shock  impingement  conditions  are  set 
with  PI =500  Pa  and  normal  shock  acting  at  the  center  of  the  gage  head.  This  way  we 
know  the  gage  will  work  under  less  stringent  but  more  realistic  impingement  conditions. 


Table  5-3.  Test  conditions  for  consideration  for  CF3. 


Time 

Tw 

Max  Acc  (g’s) 

Max  Mom 

LENS  I 

10  ms 

1000  Pa 

1 

M10NS 

LENS  II 

50  ms 

2000  Pa 

1 

M10NS 

LENS  48  in. 

5  ms 

150  Pa 

1 

M10NS 

VTSST 

12  s 

85  Pa 

0 

M3. 7  NS 

Now  that  load  conditions  have  been  established,  we  can  begin  to  narrow  the 
design  space.  We  start  with  the  geometric  constraints  on  the  design  variables.  From  the 
previous  gages  we  see  that  the  three  design  variables  are  beam  length,  thickness,  and 
width.  By  choosing  the  FC-8  bellows  and  our  full  bridge  gage  configuration,  we  have 
effectively  frozen  the  design  width  to  0.12  in.  This  has  eliminated  a  design  variable.  The 
two  remaining  variables  are  constrained  by  the  FC8  bellows,  manufacturing,  and  a 
reasonable  maximum  on  beam  length.  This  has  considerably  narrowed  our  overall  design 
space.  These  constrained  conditions  are  presented  in  Table  5-4. 


Table  5-4.  CF3  design  variable  constraints. 


Min  (m) 

Min  (in) 

Max  (in) 

Max  (m) 

Strain  gages 

0.003048 

0.12 

w 

0.12 

0.003048 

FC8  Bellows  ID 

Manufacturing 

0.000254 

0.01 

th 

0.12 

0.003048 

FC8  Bellows  ID 

FC8  +  Lee  Plug 

0.012065 

0.475 

L 

2 

0.0508 

Model  space 
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The  physical  design  space  can  now  be  shown  graphically.  Figure  5-2  shows  the 
range  of  possible  designs  in  two  ways.  First,  on  the  left,  we  see  the  four  extremes.  For 
example,  the  top  left  gage  is  the  shortest-thinnest  possible  beam,  while  the  bottom  right 
shows  the  longest-thickest  possible  beam.  This  same  space  is  expressed  in  a  different 
way  with  the  figure  on  the  right.  In  this  case,  the  largest,  thickest  beam  is  split  in  half. 
Then  the  range  of  geometrically  allowable  thickness  and  length  are  each  divided  into  25 
steps.  In  other  words,  each  of  the  625  “pixels”  represents  a  completely  different  gage 
design.  As  we  look  down  the  rows,  the  gages  get  longer.  Going  from  left  to  right,  the 
gages  get  thicker.  The  entire  discretized  space  represents  our  geometrically  constrained 
design  space.  This  concept  is  very  important  and  will  be  the  driving  picture  to 
understanding  the  total  optimization  procedure. 


Figure  5-2.  CF3  geometrically  constrained  design  space. 

Now  that  we  have  successfully  constrained  and  discretized  the  design  space,  we 
need  to  constrain  the  state  variables.  We  start  with  seven  state  variables:  con  and  s,  5, 
m  and  eg,  and  AT.  The  first  two  make  up  the  first  mode  shape  response  of  the  gage. 
Since  we  are  still  using  an  FC-8  oil-filled  gage,  we  can  reasonably  assume  that  the 
damping  coefficient  should  remain  around  0.043.  The  output  strain  and  displacement  are 
self-explanatory.  The  mass  and  center  of  gravity  of  the  gage  help  define  how  it  will  react 
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to  inertia.  Since  we  are  interested  in  keeping  the  magnitude  of  this  response  low 
compared  to  the  wall  shear  output,  we  will  calculate  an  equivalent  inertial  shear  based  on 
maximum  g-loading.  Finally,  the  AT  is  the  change  in  Celsius  for  the  top  set  of  strain 
gages.  Thus,  with  these  adjustments,  we  have  simplified  our  problem  to  five  state 
variables:  f,  s,  8,  xa,  AT.  Next  we  must  constrain  these  variables. 

We  can  start  with  the  time  constraint  on  minimum  natural  frequency  required. 
Since  our  minimum  run  time  is  10  ms,  we  must  ensure  that  the  gage  will  settle  within  that 
time.  Figure  5-3  shows  that  a  1000  Hz  gage  will  meet  that  requirement  by  settling  within 
8.5  ms.  Thus  we  can  safely  design  a  gage  with  a  minimum  1000  Hz  natural  frequency  to 
overcome  transient  effects.  There  is  no  maximum  on  natural  frequency,  so  we  move  to 
the  next  state  variable. 


Figure  5-3.  8.5  ms  settling  time  for  1000  Hz  gage  for  10  ms  test  run. 

The  next  state  variable  for  consideration,  s,  must  be  examined  against  each  of  the 
three  different  load  cases  of  interest.  The  gage  must  produce  some  minimal  strain  output 
even  under  the  minimum  wall  shear  of  1000  Pa.  Therefore,  we  set  that  minimum  to  0.5 
ps  for  this  condition.  This  is  very  low,  but  provides  a  reasonable  start  point.  All  of  the 
analysis  of  the  gage  assumes  that  it  operates  within  the  proportional  limit  of  the  material. 
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For  INVAR,  that  limit  is  about  1.18E11  Pa.  Therefore,  we  cannot  exceed  this  limit  even 
during  overshoot  transient  conditions  for  our  worst  load  conditions,  Case  2  and  Case  3. 
Since  the  overshoot  is  just  under  twice  the  step  response,  we  can  safely  check  both 
conditions  using  twice  the  expected  step  load. 

The  displacement,  5,  has  no  minimum.  However,  to  keep  the  gap  size  to  a 
minimum,  we  set  this  distance  to  a  maximum  of  0.005  in.  To  ensure  that  the  head  does 
not  touch  the  housing  even  during  the  overshoot  conditions  of  Cases  2  and  3,  we  must 
check  the  displacement  under  twice  these  load  conditions. 

Since  each  pixel  represents  a  new  gage  design,  there  are  625  different  total  masses 
and  center  of  gravity  locations.  These  properties  will  be  converted  to  calculate  an 
equivalent  wall  shear  based  on  a  1  g  inertial  loading.  This  technique  was  presented  in 
gage  characterization.  This  equivalent  acceleration  shear  magnitude  can  be  compared  to 
our  xw  min  =  1000  Pa.  Since  we  have  greatly  increased  the  overall  size  and  density  of  the 
gage,  we  set  the  maximum  xa  to  25%  of  xwmin  or  250  Pa. 

The  final  state  variable  for  consideration  is  the  temperature  variation  at  the  top  set 
of  strain  gages,  AT.  Based  on  the  resulting  output  from  the  VT  testing,  we  decided  that 
we  wanted  to  keep  the  temperature  change  at  this  set  of  strain  gages  below  1  deg  C.  This 
particular  variable  requires  longer  FEM  run  times  for  analysis,  so  will  be  examined  for 
the  final  design  selected.  Intuitively  by  making  the  beam  length  longer  and  thicker,  more 
diffusion  can  be  achieved.  This  will  play  in  choosing  the  final  design  dimensions.  The 
constrained  state  variables  are  presented  in  Table  5-5.  Now  we  can  begin  the  careful 
process  of  examining  how  each  of  the  625  gage  designs  responds  to  each  of  our  load 
conditions. 

The  optimization  process  uses  four  separate  sweeps  of  all  625  designs.  First  the 
axisymmetric  bellows  model  is  used  to  determine  first  mode  natural  frequency.  Second, 
a  3-D  model  is  used  to  determine  the  maximum  sy  and  maximum  8X  due  to  a  2000  Pa 
load.  Third,  a  3-D  model  is  used  to  determine  the  maximum  sy  and  maximum  8X  due  to 
the  worst  case  moment  load  condition  of  a  Mach  10  normal  shock  impinging  on  the 
center  of  the  gage  head.  Notice  that  the  load  conditions  for  the  wall  shear  and  moment 
case  are  extreme.  The  beauty  of  ensuring  that  the  output  strains  are  below  the 
proportional  limit  is  linearity.  In  other  words,  once  we  have  a  strain  output  for  the  2000 
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Pa  case,  as  long  as  we  stay  within  the  proportional  limit,  the  strain  output  for  the  1000  Pa 
case  is  simply  half.  Finally,  an  equivalent  inertial  magnitude  sweep  is  conducted. 


Table  5-5.  CF3  constrained  state  variables. 


LOAD 

min 

SV 

max 

LOAD 

0.85  Min  t 

1000  Hz 

f 

NA 

Min  tw 

0.5  pe 

8 

8  pi 

2  tw  max 

0 

8 

8  pi 

2  Mo  max 

0 

Ta 

0-25  twmin 

Max  acc 

0 

8 

5  mil  (gap) 

2  xw  max 

0 

8 

5  mil  (gap) 

2  Mo  max 

0 

AT 

1  deg  C 

qw(t) + 1 

The  first  state  variable  considered  is  first  mode  natural  frequency,  f.  Figure  5-4 
shows  how  this  state  variable  varies  with  the  design  space.  The  left  figure  shows 
unbound  natural  frequency  in  Hertz.  The  design  space  shows  that  the  beam  length 
increases  from  top  to  bottom,  and  beam  thickness  increases  from  left  to  right.  We  can  tell 
as  we  choose  designs  that  go  up  and  right,  or  get  shorter  and  thicker,  we  can  achieve 
higher  frequencies.  In  order  to  bind  this  design  space,  the  figure  to  the  right  shows  the 
effect  of  eliminating  any  gages  that  have  frequencies  less  than  1000  Hz.  The  designs  that 
do  not  meet  this  minimum  requirement  are  removed  as  possible  designs.  Clearly  the 
limit  on  this  state  variable  has  greatly  reduced  the  number  of  possible  designs  and  has 
effectively  constrained  us  to  a  shorter  thicker  gage  than  would  be  demanded  just  due  to 
geometric  constraints. 
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Figure  5-4.  CF3  first  mode  natural  frequency. 

The  next  state  variable  examined  is  strain.  Figure  5-5  shows  how  each  design 
responds  to  the  minimum  expected  wall  shear,  1000  Pa.  The  unbound  contour  values  of 
the  left-hand  figure  propose  ridiculous  strain  outputs  where  the  material  would  really 
have  failed.  These  impossible  designs  are  eliminated  in  the  right  hand  figure  by  applying 
the  proportional  limit  criterion  for  an  overshoot  2000  Pa  load,  i.e.  4000  Pa.  This 
effectively  eliminates  extremely  long  and  thin  beam  designs.  Of  course  these  were 
already  eliminated  by  the  minimum  frequency  requirement.  However,  by  maintaining 
careful  bookkeeping  of  each  state  variable  we  can  confidently  ensure  we  understand  the 
gage  output  within  these  extreme  load  conditions.  Additionally,  the  right  hand  figure 
shows  the  impact  of  eliminating  low  strain  designs.  In  direct  opposition  to  frequency 
requirements,  this  minimum  strain  has  pushed  the  design  space  away  from  the  shortest 
thickest  available  designs.  The  minimum  strain  requirement  actually  pushes  the  design 
space  down  and  to  the  left,  i.e.  we  desire  longer  thinner  beams.  This  direct  contrast  to 
frequency  considerations  was  introduced  in  the  design  of  CF1  and  continues  to  be  a 
driver  for  designing  direct  measure  gages. 
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Figure  5-5.  CF3  strain  output  due  to  1000  Pa  wall  shear. 


Displacement  must  be  verified  to  ensure  the  intended  5  mil  gap  will  suffice. 
Since  the  second  optimization  sweep  allows  us  to  predict  the  maximum  8X  for  each 
design  we  can  check  the  effects  of  geometry  on  this  state  variable.  Figure  5-6  shows  the 
expected  displacements  for  our  2000  Pa  case.  In  order  to  eliminate  designs  where 
overshoot  loading  would  allow  the  head  to  hit  the  housing,  we  again  remove 
unacceptable  designs  in  the  right  hand  figure.  As  with  the  strain  proportional  limit,  the 
minimum  natural  frequency  limit  has  already  eliminated  these  potential  designs. 
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Figure  5-6.  CF3  displacement  due  to  2000  Pa  load. 


Two  of  the  previous  state  variables,  e  and  5,  must  be  checked  against  the  moment 
loading  condition.  This  moment  sweep  did  not  change  the  already  constrained  design 
space.  Therefore,  we  will  save  the  results  from  this  sweep  for  further  discussion  after 
looking  at  inertial  effects. 

The  final  design  sweep  is  examines  the  relative  effect  of  inertia  to  wall  shear. 

This  is  done  by  calculating  an  equivalent  inertial  shear  due  to  1  g  acceleration.  Figure 
4-4  shows  this  process.  The  effect  of  even  a  small  1  g  load  is  pronounced  with  these 
INVAR  designs  due  to  the  increased  mass.  Figure  5-7  shows  that  all  of  the  625  designs 
fall  between  20-30%  of  the  expected  1000  Pa  wall  shear.  In  order  to  keep  our  design 
within  some  reasonable  range,  the  design  was  limited  to  provide  at  most  25%  equivalent 
wall  shear  due  to  1  g  acceleration.  Looking  at  the  right  hand  image  in  Figure  5-7  we  see 
the  impact  on  the  design  space  due  to  this  limitation.  Inertial  requirements  force  the 
design  space  up  and  to  the  left.  This  is  intuitively  obvious  since  we  want  a  thinner 
shorter  beam  in  terms  of  this  state  variable. 
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Figure  5-7.  CF3  equivalent  acceleration  shear. 


Now  we  can  combine  the  total  effect  of  all  design  and  state  variable  constraints. 
Figure  5-8  shows  each  of  the  previously  considered  state  variables  and  the  contribution  to 
limiting  design  space.  Additionally,  the  effect  of  moment  on  strain  and  displacement  is 
provided.  From  this  we  see  that  minimum  natural  frequency  requirement,  minimum 
strain  output,  and  maximum  inertia  dominate  for  the  LENS  design  cases.  The 
culmination  of  this  effort  is  to  narrow  the  design  to  one  that  produces  the  most  strain 
output  with  the  smallest  temperature  gradient. 

The  final  design  and  state  variable  constrained  design  space  is  shown  in  Figure 
5-9.  This  is  significantly  smaller  than  the  original  space  permitted  due  simply  to 
geometrical  constraints.  This  figure  illustrates  the  impact  of  the  major  state  variables  on 
design.  First,  the  natural  frequency  requirements  push  the  design  up  and  to  the  right. 
Second,  minimum  strain  requirements  counteract  this  by  pushing  down  and  left.  The 
trade-off  between  these  variables  was  already  well  understood.  The  addition  of  inertial 
impact  and  thermal  considerations  is  interesting.  The  desire  to  make  the  beam  longer  and 
thicker  to  allow  greater  thermal  diffusion  is  directly  opposed  by  the  desire  to  make  the 
beam  lighter.  In  total,  these  four  design  variables  have  contracted  the  design  space  to  a 
small  band  of  designs  that  will  meet  all  of  the  design  requirements  imposed  by  LENS  test 
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conditions.  Since  we  would  like  to  achieve  the  most  strain  possible  with  the  smallest 
thermal  gradients,  the  circled  pixel  represents  the  chosen  configuration  for  the  CF3  gage. 
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Figure  5-8.  CF3  state  variable  constraints  on  design  space. 


Figure  5-9.  CF3  fully  constrained  design  space  with  expected  strain  output. 
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5.1.2  CF4  Gage  Design 

After  completing  the  design  for  the  INVAR  gage,  we  recognized  the  desire  to 
reduce  weight.  The  true  driver  for  using  INVAR  was  to  avoid  thermal  problems 
associated  with  the  long  transient  heat  flux  effects  in  the  Virginia  Tech  supersonic  tunnel. 
Therefore,  we  decided  to  make  an  aluminum  gage  that  would  meet  all  of  the  LENS 
requirements.  Fortunately,  the  design  optimization  tools  developed  for  CF3  could  be 
readily  employed  to  develop  an  aluminum  gage.  The  same  geometric  and  state  variable 
constraints  apply.  These  were  presented  in  Table  5-4  and  Table  5-5  respectively.  The 
same  four  sweeps  were  conducted:  1)  natural  frequency  2)  wall  shear  impact  on  s  and  5 
3)  moment  impact  on  s  and  8  and  4)  equivalent  acceleration  shear.  Using  this  process  we 
can  look  directly  at  the  effect  of  material  choice  on  design  space.  Additionally,  although 
all  state  variables  are  checked,  we  can  limit  the  presentation  of  optimization  to  the  four 
design  space  drivers  shown  in  Figure  5-9. 

Starting  with  the  same  minimum  requirement  for  natural  frequency,  we  can 
compare  the  aluminum  design  space  to  that  of  INVAR.  This  is  shown  in  Figure  5-10. 

The  major  trend  to  push  up  and  to  the  right  remains  the  same.  Additionally,  the 
frequencies  are  in  similar  range  for  both  materials.  Also,  there  are  more  allowable 
configurations  for  aluminum. 

The  next  state  variable  is  the  strain  output.  The  greater  ductility  of  aluminum 
significantly  impacts  the  design  here.  Figure  5-1 1  shows  that  even  with  10  times  the 
required  minimum  strain  output,  the  design  space  is  still  bigger  than  that  of  INVAR.  We 
can  expect  much  larger  strain  outputs  from  an  optimized  aluminum  skin  friction  gage 
than  from  CF3.  This  is  an  additional  advantage  to  limiting  apparent  strain  due  to  inertia. 
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Figure  5-10.  CF4  frequency  impact  on  design  space,  aluminum  vs.  invar. 
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Figure  5-11.  CF4  strain  response  to  wall  shear. 


The  last  state  variable  for  comparison  is  equivalent  acceleration  shear.  This  is  the 
reason  we  switched  from  INVAR  back  to  aluminum.  Clearly  the  change  is  beneficial. 
Now  the  range  of  expected  equivalent  acceleration  shear  has  decreased  to  15-18%  of 
1000  Pa  for  a  1  g  acceleration.  In  fact  there  is  no  limitation  on  design  space  based  on  our 
previous  criteria  for  CF3. 
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Figure  5-12.  CF4  equivalent  acceleration  shear  as  percentage  of  wall  shear. 


Now  we  can  combine  the  state  variable  impact  on  design  space  as  previously 
performed  for  CF3.  This  is  compared  to  the  INVAR  design  space  in  Figure  5-13.  Notice 
how  many  more  possible  solutions  there  are  using  aluminum.  However,  we  must  be 
careful  to  verify  low  thermal  effects.  This  was  accomplished  in  Chapter  4:  Gage 
Characterization.  At  this  point,  we  use  engineering  intuition  to  choose  our  design  for 
CF4.  Since  the  inertial  effects  are  low,  we  can  choose  the  thickest,  longest  available 
design.  Aside  from  the  expected  thermal  benefit,  this  skin  friction  gage  could  be  outfitted 
for  full  3-D  measurements.  In  other  words,  we  could  place  strain  gages  on  all  beam  faces 
and  measure  both  axial  and  spanwise  shear  with  this  design. 
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Figure  5-13.  CF4  constrained  design  space. 


5.2  Optimization  for  Blow  Down  Facility  Testing 

Since  our  test  verification  facility  at  Virginia  Tech  is  quite  different  from  LENS, 
it  is  interesting  to  see  what  happens  to  the  design  space  with  these  different  constraints. 
Since  we  know  that  longer  test  times  have  a  detrimental  impact  on  temperature  gradients, 
we  will  apply  an  INVAR  design.  Rather  than  step  through  each  output  variable,  the 
entire  design  solution  is  presented  in  Figure  5-14.  The  design  conditions  have  a  dramatic 
impact  on  which  parameters  are  drivers.  The  short  duration  facilities  settling 
requirements  forced  the  design  space  up  and  right.  Here  the  low  wall  shear  has  forced 
the  design  space  down  and  to  the  left.  The  other  main  driver  for  the  blow  down  facility 
comes  from  the  desire  to  avoid  railing  during  moment  testing.  The  low  acceleration 
experienced  in  the  blow-down  facilities  removes  inertia  as  a  constraining  variable. 
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Figure  5-14.  Blow  down  facility  constraints  on  design  space. 


Now  that  we  have  seen  how  the  design  space  has  been  affected  by  each  state 
variable,  we  can  combine  the  effects  and  look  at  the  expected  strain  output  for  design 
space  under  an  85  Pa  wall  shear.  Figure  5-15  shows  that  the  design  can  be  considerably 
longer  and  thinner  than  CF3.  In  fact,  because  of  the  less  stringent  requirement  on 
vibrational  frequency,  even  longer  designs  could  be  considered. 
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Figure  5-15.  Blow  down  facility  expected  wall  shear  strain  output. 
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5.3  Optimization  for  a  Typical  Flight  Test 


The  next  test  condition  for  consideration  is  a  flight  test.  As  with  the  blow  down 
facility,  the  long  test  time  makes  INVAR  an  attractive  material  choice.  Figure  5-16 
shows  the  overall  impact  of  constraints  on  the  state  variables.  As  with  the  blow  down 
facility,  this  test  condition  does  not  limit  the  design  due  to  vibrational  frequency.  In  this 
case  however,  the  higher  expected  acceleration  returns  as  a  limiting  factor.  Figure  5-17 
shows  the  total  effect  of  all  constraints  on  design  space  and  the  expected  wall  shear 
induced  strain  output. 


Figure  5-16  Flight  test  constraints  on  design  space. 


Figure  5-17  Flight  condition  expected  wall  shear  strain  output. 
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5.4  Optimization  Summary 


The  optimization  technique  developed  through  this  chapter  is  a  straightforward 
parametric  analysis.  First,  the  material  choice  was  made  and  the  design  variables  were 
limited  to  a  constrained  length  and  thickness.  Then,  625  different  gages  within  that 
design  space  were  virtually  built  and  tested.  Once  we  had  the  gage  performance 
characteristics  for  every  gage,  we  compared  their  performance  to  test  facility 
specifications.  Finally,  a  design  that  fit  within  those  constraints  and  provided  the  desired 
performance  was  selected. 

Figure  5-18  shows  the  inputs  and  corresponding  state  variables  of  interest  as  well 
as  generalized  impact  on  design  space.  Notice  that  the  modeling  techniques  chosen  in 
Chapter  4,  Gage  Characterization,  have  been  applied  625  times  for  INVAR  then  another 
625  times  for  aluminum.  The  only  difference  in  making  the  design  decision  was  how  the 
gage  met  the  performance  constraints  of  each  facility.  The  heat  flux  model  is  only 
applied  after  gage  selection  due  to  the  long  computational  time  required  for  this 
investigation. 


Figure  5-18.  Optimization  overview. 
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6  Simulations 


The  purpose  of  developing  models  was  twofold:  design  optimization  and  virtual 
testing.  A  full  simulation  connects  the  flow,  gage,  and  instrumentation  models. 
Simulations  can  be  either  experimental  or  computation,  and  both  approaches  have  been 
used  in  this  work.  The  experiments  were  discussed  in  Chapter  3.  Through  simulation  we 
can  answer  important  questions.  For  example,  the  simulation  developed  in  the  Chapter  2 
Theoretical  Approach  answered  how  to  configure  the  Wheatstone  bridges  to  separately 
measure  wall  shear  and  moment  while  being  insensitive  to  other  loads. 

Several  important  questions  required  answering  throughout  this  research.  These 
were  attacked  through  simulation.  Each  of  the  following  sections  poses  an  important 
physical  question  and  the  simulation  solution.  Unless  noted  the  instrumentation  is 
assumed  ideal  so  that  the  gage  output  is  exactly  what  is  measured. 

6.1  Instrumentation 

The  sources  of  apparent  strain  for  this  type  of  gage  were  discussed  in  Chapter  4: 
Gage  Characterization.  Additionally  we  must  investigate  other  possible  causes  of  error. 
Instrumentation  provides  the  other  most  likely  reason  for  measurement  error.  These 
sources  come  from  strain  gage  manufacture  and  alignment,  lead  wire  resistance,  and 
amplifier  resistances.  Table  6-1  lists  these  error  sources  and  maximum  impact  which  will 
be  applied  to  determine  overall  system  error.  All  of  these  errors  are  quoted  by  the 
manufacturer  except  strain  gage  misalignment.  These  are  conservatively  high  estimates 
based  on  FEM  strain  distribution. 


Table  6-1.  Instrumentation  error  sources. 37 


Maximum  error 

Strain  gage  Rg 

+/-  0.2  %  =  0.7  Q 

Strain  gage  vertical  misalignment  1  mm  offset 

5%  of  total  ps 

Strain  gage  angular  misalignment  5  deg 

1%  fl£ 

Lead  wire  resistance  3  ft  0.35  Q/m 

0.32  Q 

Meter  resistance  5  %  of  100  kQ 

5  kQ 

Power  supply  resistance 

0.1  Q 
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As  a  worst  case  scenario,  all  of  the  strain  gage  errors  were  applied  to  one  side  of 
the  half  bridge  to  determine  the  maximum  error  impact  on  total  voltage  output. 
Accounting  for  this  error  and  lead  wire  resistance  we  can  find  the  expected  real  output 
voltage  using  Equation  6-1.  Using  the  nominal  output  conditions  for  CF3  from  Chapter 
4,  yields  a  5.3  %  variation  of  output  voltage  at  the  root. 

i?i  (. r2  +  rl  )—  r2  X  +  rl  )  - 

AV  =  Vm  Gain  - — - — - - - -  Equation  6-1 38 

X  +R2+  2Rl  )  +  (/?,  +R2+2Rl  XA Rl  +  A R2 ) 


Now  accounting  for  the  additional  error  of  impedance  from  power  supply  and 
meter,  we  can  determine  the  total  error  output  for  the  system.  This  is  accomplished  by 
multiplying  Equation  6-1  by  the  following  factor  where  Rs  is  supply  resistance,  Rm  is 
meter  resistance,  and  RB  is  bridge  resistance,  calculated  to  be  350.175  Q.  The  resulting 
total  error  on  output  voltage  due  to  these  maximum  error  conditions  is  5.6  %  at  the  root 
for  CF3. 


factor  = 


1 

i 

R 

R„ 

1  +  -5- 

1  +  -^- 

L 

Rm  _ 

Equation  6-2 


If  there  is  no  strain  gage  vertical  or  angular  misalignment,  the  expected  error  is 
1.7%  from  lead  wire  resistance  alone  and  a  total  instrumentation  error  of  2.0%.  This 
analysis  has  shown  that  the  maximum  expected  instrumentation  errors  are  relatively  low. 
Most  of  the  error  comes  from  the  high  estimate  of  alignment  error.  Other  instrumentation 
error  is  due  to  lead  wire  resistance,  which  is  limited  by  keeping  the  wires  short.  In  the 
sections  that  follow,  instrumentation  error  is  assumed  zero  unless  noted. 


6.2  Dynamic  Response 

The  entire  measurement  process  for  these  types  of  gages  is  dynamic.  Even  an 
ideal  step  load  produces  transient  response.  This  aspect  of  the  instrumentation  brought 
up  several  issues  that  were  resolved  through  simulation. 
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6.2.1  Minimum  Natural  Frequency 


One  of  the  first  hurdles  for  this  type  of  gage  is  to  achieve  fast  enough  settling  time 
to  provide  accurate  measurements  even  in  the  short  run  times  encountered  in  LENS.  The 
first  question:  What  is  the  minimum  natural  vibration  frequency  necessary  for  a  settling 
within  a  10  ms  test  run?  The  conservative  answer  of  1000  Hz  was  presented  in  Figure 
5-3.  But  this  is  only  a  part  of  our  total  simulation  solution. 

6.2.2  Filter  Setting 

This  brought  up  another  important  question.  What  Butterworth  filter  should  be 
used  to  provide  the  best  signal  output?  If  we  assume  a  simple  step  load  of  1000  Pa,  and 
an  exactly  calibrated  gage  with  ideal  instrumentation,  we  can  simulate  the  signal  response 
for  several  different  cut-off  frequencies  as  in  Figure  6-1.  Here  we  see  that  the  100  Hz 
filter  is  the  bare  minimum  to  try  and  measure  this  gage  reaction.  A  1000  Hz  or  10  kHz 
filter  is  preferable.  As  a  note,  wideband  or  raw  signal  response  could  be  used.  However 
it  would  produce  aliasing  if  we  wished  to  perform  signal  analysis  in  the  frequency 
domain. 


Figure  6-1.  Expected  signal  output  using  several  cut-off  frequencies. 
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6.2.3  Gage  Reaction  to  Heaviside  load 

We’ve  seen  how  the  gage  should  respond  to  a  simple  step  load.  But  what  about  a 
load  that  is  more  representative  of  test  conditions?  Additionally,  we  have  to  consider  the 
impact  of  imperfect  instrumentation.  So  we  have  to  answer  two  questions:  1)  how  does 
the  gage  response  look  given  a  realistic  load?  and  2)  what  kind  of  error  can  we  expect 
with  real  instrumentation  error? 

A  full  simulation  was  developed.  Figure  6-2  shows  a  simplified  Heaviside 
forcing  function  approximation  to  the  time  varying  wall  shear  measured  in  the  blow¬ 
down  facility  at  Virginia  Tech.  The  beam  model  for  displacement  and  strain  is  a 
cantilever.  The  voltage  output  comes  from  a  realistic  instrumentation  system  that 
accounts  for  strain  gage  vertical  and  angular  misalignment,  lead  wire  resistance,  and 
power  supply  impedance. 


SST  Shear:  ANSYS  iw  Load  (Plate) 


Figure  6-2.  Heaviside  approximation  to  Virginia  Tech  time  varying  wall  shear. 

A  screen  shot  near  the  end  of  the  steady  load  of  this  simulation  is  presented  as 
Figure  6-3.  In  the  top  left  image,  we  see  the  wall  shear  deck.  In  the  top  right  is  the  gage 
displacement  with  strain  gages  slightly  offset  near  the  head.  The  bottom  left  shows  the 
strain  measurements  at  each  strain  gage  location.  Finally,  the  bottom  right  shows  the 
resulting  expected  voltage  output  for  this  load  condition. 


91 


In  answer  to  our  first  question,  the  gage  response  to  this  single  load  input  was  as 
we  expected.  Notice  that  no  filter  was  applied.  This  is  evidenced  by  the  rapid  overshoot 
captured  in  the  bottom  left  figure.  So  for  long  run  times  as  in  the  VT  supersonic  tunnel, 
we  can  overcome  the  transient  response  by  application  of  a  low-pass  filter. 

In  answer  to  the  second  question,  the  packaged  instrumentation  showed 
considerable  robustness.  The  simulation  showed  that  despite  the  introduction  of  multiple, 
but  reasonable,  error  sources  within  the  instrumentation,  the  overall  system  error 
excepting  strain  gage  alignment  was  around  1%.  The  simulation  also  showed  that  even 
with  slightly  exaggerated  misalignments  of  the  strain  gages,  the  measurements  were 
within  5%  of  actual.  This  was  good  news  since,  such  a  large  error  in  manufacturing 
would  have  had  the  gages  offset  by  1/2  gage  length.  Fortunately,  none  of  our  gages  are 
offset  to  the  extent  that  could  be  detected  by  the  naked  eye,  which  is  much  less  than 
1/1 0th  of  a  gage  length. 


Tauw  at  t=  19.8s 


IDO 


re'  50 

Cl 


CT5 


0 


-50 


i= 

05 

i _ 

"tO 

£ 


10  20 
t(s) 

Strain 


t(s) 


Displacement  vs  x  at  t=  19.8s 


30 


x  -iq'5  Vout  due  to  S 


> 


t(s) 


Figure  6-3.  Simulation  of  Virginia  Tech  pure  wall  shear  case. 
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6.2.4  Base  Acceleration  Response39 

The  loads  and  initial  conditions  evaluated  to  this  point  were  due  to  single  input 
systems.  Accelerometer  measurements  in  the  LENS  facility  showed  that  a  multiple  input 
simulation  was  necessary.  In  this  case,  we  developed  transfer  functions  to  capture  base 
acceleration  and  wall  shear  induced  motion.  The  “base”  in  this  case  is  the  physical  test 
specimen,  for  example  a  model  of  the  NASA  Aerospace  plane.  The  plane  would 
experience  acceleration  during  the  tunnel  run  similar  to  that  shown  in  Figure  6-4.  Using 
the  trapezoidal  rule,  the  velocity,  and  displacement  of  the  base  are  also  shown.  Here  we 
see  that  in  just  25  ms,  the  base  model  has  displaced  over  half  a  millimeter.  This  would  in 
turn  affect  the  skin  friction  gage.  Figure  6-5  shows  a  simple  uni-axial  block  model  of  the 
skin  friction  gage  in  this  condition. 


Figure  6-4.  Acceleration,  velocity,  and  displacement  of  base  in  LENS. 
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Figure  6-5. 1-D  kinetics  diagram. 
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This  model  can  be  broken  down  into  two  linearly  independent  systems  where  the 
total  motion,  xL,  is  composed  of  the  superposition  of  motion  due  to  base  acceleration,  xu? 
and  the  motion  due  to  shear  load,  xs.  This  is  Equation  6-3.  The  first  component  of 
displacement  considered  is  due  to  a  step  wall  shear,  S(t).  Equation  6-4  shows  the 
corresponding  ordinary  differential  equation  (ODE).  The  coefficients  Du,  Kn,  and  mT 
are  gage  properties  and  L  is  length  along  the  beam.  The  next  component  of  displacement 
is  due  to  the  base  motion  shown  in  Figure  6-4.  The  ODE  for  this  component  is  Equation 
6-5. 

XL  =  ,VS  +  xu  Equation  6-3 

mTxs  +  Duxs  +  Kuxs  =  S(t )  Equation  6-4 

mTXu  +  Dnxu  +  KnXu  =  Dnu  +  Knu  Equation  6-5 

We  are  interested  in  the  relative  motion  of  the  gage  from  the  base  reference 
frame.  Therefore,  we  have  to  recast  Equation  6-5.  Consider  that  the  relative  motion  of 
the  gage  due  to  the  base  motion  is  z  =  xu  -  u,  where  xu  is  the  total  gage  displacement  due 
to  base  motion  and  u  is  the  base  displacement.  Substituting  this  relative  motion  equation 
into  Equation  6-5  and  solving  yields  Equation  6-6. 
mTz  +  Dnz  +  Knz  =  -mTii  Equation  6-6 

The  corresponding  transfer  functions  for  Equations  6-4  and  6-6  are  Equations  6-7 
and  6-8  respectively.  Using  the  Matlab  ®  control  system  toolbox  for  linear  time 
invariant  systems,  we  can  simulate  the  system  response  to  the  wall  shear  and  acceleration 
inputs.  Figure  6-6  shows  the  result  of  this  simulation.  In  the  top  left  figure,  the  response 
to  a  pure  wall  shear  of  2000  Pa  acting  on  the  CF3  gage  is  modeled.  The  predicted 
displacement  of  about  4.5  pm  is  of  the  right  order  of  magnitude  of  the  FEM  predicted  3.1 
pm.  The  top  right  figure  shows  base  motion  induced  displacement.  In  this  case,  the 
dynamic  model  predicts  a  maximum  absolute  displacement  of  about  0.3  pm, 
corresponding  to  lg  acceleration.  This  matches  the  combined  solid  mechanics/  FEM 
prediction  of  0.3  pm. 
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Equation  6-7 
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It  is  important  to  note  that  this  1-D  dynamic  response  model  is  not  intended  to 
accurately  predict  displacement.  The  1-D  simplification  causes  this  simulation  to  slightly 
overpredict  the  displacement  due  to  wall  shear.  This  simulation  was  intended  to  show  a 
relative  dynamic  system  response  to  two  separate  inputs.  This  is  shown  in  the  bottom 
image  of  Figure  6-6.  This  simulation  of  the  total  signal  output  provides  several  insights. 
First,  the  total  response  is  not  static.  The  real  measurement  signal  is  going  to  have 
expected  oscillations.  Second,  if  we  look  just  at  the  average  signal  values  over  the  test 
window,  the  pure  shear  output  is  not  too  different  from  the  total  signal  output.  In  fact,  the 
variation  is  just  1.2  %  and  that’s  with  this  simplified  model. 


Figure  6-6.  Simulated  1-D  gage  response  to  multiple  inputs. 

Thanks  to  this  investigation,  acceleration  considerations  could  be  justifiably 
simplified.  Rather  than  trying  to  track  the  overall  effect  on  the  signal,  we  scaled  back  to 
looking  at  expected  order  of  magnitude  effects.  Thus  we  developed  the  equivalent 
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acceleration  shear.  This  quick  approach  allows  a  convenient  method  to  design  gages 
constrained  by  designer  imposed  inertial  limits.  Additionally,  the  simple  solid  mechanics 
model  can  more  accurately  predict  the  relative  strain  output. 

6.3  Thermal  response 

The  thermal  issue  uncovered  during  the  Virginia  Tech  Shear  test  spawned  several 
questions.  Was  there  a  direct  relation  between  temperature  and  voltage  output?  This  was 
eliminated  as  a  possibility  through  calibration  testing.  This  meant  there  was  either 
something  wrong  with  the  instrumentation,  thermal  gradients  were  causing  a  problem,  or 
the  temperature  distribution  varied  through  the  thickness  of  the  beam.  Each  of  these 
possibilities  was  examined  carefully  through  simulation.  The  simulations  started  from  a 
simple  1-D  finite  difference  model  (FDM),  to  2-D  FDM,  to  Axisymmetric  FEM,  and 
finally  a  full  3-D  FEM.  Additionally,  the  boundary  conditions  increased  in  complexity 
from  a  varying  wall  temperature  to  a  time  and  spatially  varying  heat  flux. 

Starting  with  the  1-D  case,  we  can  examine  expected  thermal  changes  due  to  step 
loads.  The  method  of  verification  is  the  semi-infinite  body  Gaussian  error  function 
distribution.  Two  types  of  analytical  solutions  are  available.  One  is  for  a  step 
temperature  condition.  The  second  is  for  a  step  heat  flux.  Both  only  consider  1  -D  axial 
changes,  so  for  the  purpose  of  this  model  “x”  is  used  to  describe  this  direction  and 
denotes  distance  from  the  surface.  This  direction  is  negative  y  in  our  Cartesian  flow 
frame. 

6.3.1  Step  Temperature  Boundary  Condition 

The  step  boundary  condition  provides  the  basis  for  determining  heat  flux  from 
static  wall  temperature.  We  verified  our  1st  order  time  2nd  order  space  forward  time 
central  space  (FTCS)  solution  to  the  heat  flux  equation.  An  aluminum  semi-half  body, 
a=8.6E-5  m  /s,  initially  at  290K  has  a  step  surface  temperature  change  to  400  K.  This 
step  boundary  condition,  TWO=290  K  and  Tw2=400K,  is  applied  for  9.0  s.  Table  6-2 
shows  the  governing  equations  for  the  model  and  analytic  solutions.  For  the  FTCS,  n  is 
time  index,  i,  is  the  spatial  index.  The  FTCS  model  met  the  stability  criteria  of  d  fA3.5. 
The  model  results  were  well  matched.  The  final  result  at  9.0s  is  shown  in  Figure  6-7. 
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Table  6-2.  Thermal  step  equations. 
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Figure  6-7.  Gaussian  analytic  and  1-D  FTCS  thermal  step  solution  at  9.0  s. 


6.3.2  Step  Heat  Flux  Boundary  Condition 

The  step  heat  flux  boundary  condition  is  useful  especially  for  a  first  estimate  of 
thermal  effects.  For  example,  we  have  an  estimate  of  the  constant  heat  flux  value  in  the 
LENS  facility.  This  verification  uses  an  aluminum  semi-body  initially  at  287  K  with  a 
sudden  heat  flux  of  100  kW/m  .  This  condition  is  applied  for  0.9s.  The  equation  for  this 
condition  is  presented  in  Table  6-3. 
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Table  6-3.  Analytic  solution  for  sudden  heat  flux. 


Analytic42 
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A  comparison  is  made  between  the  FTCS  model  and  the  analytic  model.  The 
temperature  distributions  for  both  solutions  after  0.9  s  are  shown  in  Figure  6-8.  As  with 
the  thermal  step  solution,  the  numerical  approximation  lies  directly  on  top  of  the  analytic 
solution. 


Figure  6-8.  Gaussian  analytic  and  1-D  FTCS  heat  flux  step  solution  at  9.0  s. 

Heat  flux  determines  the  temperature  distribution.  Since  we  only  have  a 
thermocouple  measurement  at  a  point  on  the  surface,  we  need  to  use  a  method  to 
approximate  heat  flux.  Several  techniques  are  available.43  Unfortunately,  none  of  these 
techniques  were  satisfactory.  A  new  method  was  developed  that  emphasized  the  test  run 
physics.  We  know  the  tunnel  sees  a  1)  starting  shock  followed  by  2)  a  transient  state  and 
finally  3)  the  steady  state.  So  rather  than  trying  to  develop  an  exact  time  varying  heat 
flux,  we  only  seek  three  load  step  conditions  that  correspond  to  these  three  states. 
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The  result  of  this  technique  is  compared  to  the  classic  Cook-Felderman 
integration  method  in  Figure  6-9.  Clearly  the  new  technique  provides  a  “clean” 
approximation  to  the  heat  flux  that  is  relatively  simple  to  input  as  a  surface  flux  in  an 
ANSYS  model. 


Figure  6-9.  Heat  flux  development  for  Virginia  Tech  supersonic  tunnel  using  thermocouple  data. 

The  validation  of  this  simplified  heat  flux  step  load  uses  an  Axisymmetric  FEM 
where  a  node  at  the  thermocouple  location  tracks  temperature.  The  result  is  shown  in 
Figure  6-10.  Even  though  the  temperature  output  is  slightly  different  than  measured,  the 
difference  is  negligible  and  worth  the  simplification  of  the  load  step  approximation. 
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Figure  6-10.  Validation  of  Axisymmetric  FEM  and  approximated  heat  flux. 

The  combined  approximate  heat  flux  and  wall  shear  are  shown  in  Figure  6-11.  It 
is  important  to  notice  that  even  these  approximations  require  10  separate  load  steps. 
However  these  load  conditions  provide  a  reasonable  approximation  to  allow  for 
determining  gage  response. 
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Figure  6-11.  Load  deck  for  Virginia  Tech  supersonic  tunnel. 


6.3.3  Thermal  Distribution  2-D  FDM 

A  1-D  model  cannot  capture  the  effects  of  thermal  diffusion  due  to  oil.  What 
does  the  thermal  distribution  look  like  using  a  2-D  FDM?  In  this  case  we  used  the 
alternating  direction  implicit  method.  Again  the  measured  static  temperature  was  set  as 
the  surface  boundary  condition.  Figure  6-12  shows  the  2-D  FDM  results  late  in  the 
supersonic  shear  run.  On  the  left  we  see  the  measured  thermocouple  data.  The  figure  on 
the  right  shows  a  slice  of  the  gage.  The  bellows  is  prominently  shown.  The  air  cavity  is 
on  the  left  and  the  oil  cavity  is  on  the  right. 
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Figure  6-12.  2-D  FDM  temperature  distribution. 

Despite  this  being  only  a  2-D  model  using  a  simple  boundary  condition,  it 
provides  important  visual  benefit.  Notice  how  the  temperature  distribution  has  diffused 
quickly  through  the  air  and  is  slowed  by  the  oil.  This  shows  that  the  oil  provides  some 
thermal  protection  at  least  from  the  air.  Additionally,  this  figure  shows  how  the 
temperature  gradients  near  the  head  strain  gage  are  much  greater  than  at  the  root.  This 
simulation  has  helped  to  provide  some  physical  insight  but  raises  other  important 
questions. 

6.3.4  2D  FDM  Result  with  Instrumentation  Error 

In  for  temperature  to  produce  an  output  voltage,  there  has  to  be  a  difference 
between  the  temperature  measured  by  the  strain  gage  at  the  front  and  rear  of  the  beam. 
The  only  way  for  this  to  occur  in  the  2D  simulation  would  be  from  strain  gage 
misalignment.  Could  an  offset  of  the  head  strain  gages  measure  enough  strain  difference 
to  cause  the  voltage  output  shown  during  testing?  The  previous  2-D  FDM  solution  to  the 
heat  equation  was  used  to  model  the  gage  response.  The  surface  boundary  condition  was 
held  to  measured  thermocouple  temperature  for  convenience.  The  strain  gages  were 
purposely  offset  in  the  model.  The  results  are  shown  in  Figure  6-13. 
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Figure  6-13.  Thermal  investigation  of  strain  gage  offset. 


In  this  screen  shot  of  the  simulation  video,  the  top  left  image  shows  the  measured 
output  for  static  temperature.  The  top  right  image  shows  the  temperature  distribution 
with  offset  strain  gages  at  the  head.  The  bottom  left  image  shows  the  strain  output  due  to 
this  purely  thermal  load.  The  bottom  right  image  shows  the  resulting  voltage  output. 

This  simulation  showed  that  even  with  matching  front  and  rear  temperature 
distributions  on  the  beam,  an  output  voltage  could  occur  with  strain  gage  misalignment. 
The  misalignment  of  this  simulation  was  purposely  exaggerated.  This  means  the  high 
voltage  measurement  from  testing  probably  came  from  an  additional  source.  This 
provided  the  impetus  for  an  axially  varying  heat  flux  investigation. 


6.3.5  Thermal  Distribution  Axisymmetric  FEM  qw(t) 

The  FDM  showed  some  potentially  interesting  thermal  gradient  issues.  In  order 
to  understand  them  better  we  need  to  increase  the  model  fidelity.  What  happens  when  we 
incorporate  fully  3-D  effects  and  use  a  realistic  heat  flux  boundary  condition?  In  order  to 
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answer  that  question,  we  started  with  an  axisymmetric  FEM  model.  Using  heat  flux  that 
is  only  dependent  on  time  as  a  boundary  condition,  we  can  look  at  a  slice  through  the 
center  of  the  gage. 


Figure  6-14.  Axisymmetric  FEM  temperature  distribution. 

Figure  6-14  is  a  screen  shot  of  the  time  varying  temperature  load.  It  shows  that 
even  with  the  oil  available  for  diffusion,  the  temperature  gradients  along  the  beam  are 
considerable.  The  temperature  range  is  relatively  small.  The  problem  is  that  when  the 
strain  due  to  shear  at  the  head  is  less  than  1  ps,  any  apparent  strain  due  to  thermal  effect 
dominates. 

6.3.6  Thermal  Distribution  3-D  qw(x,t ) 

The  final  thermal  investigation  considers  our  time  and  space  varying  heat  flux. 
What  kind  of  variation  of  temperature  distribution  is  caused  by  spatially  varying  heat 
flux?  To  answer  this  question  requires  a  full  3-D  FEM.  The  heat  flux  loading  is  shown 
in  Figure  6-15.  Each  of  the  main  steady  loads  is  highlighted  clearly  showing  the  heat 
flux  variation  with  axial  location. 
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Figure  6-15.  Time  and  space  varying  heat  flux  for  Virginia  Tech  supersonic  tunnel  pure  wall  shear. 


Next  we  can  look  at  a  full  3-D  image  of  the  real  model  and  examine  just  the 
regions  of  interest.  Although  the  3-D  FEM  solves  for  temperature  at  every  node,  we  are 
really  only  interested  in  the  temperatures  that  effect  our  strain  gages.  In  order  to  better 
see  the  thermal  distribution  at  these  important  locations,  we  look  at  planar  slices  at  the 
front  and  rear  of  the  beam.  This  concept  is  displayed  in  Figure  6-16. 
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Figure  6-16.  Planar  slices  of  3-D  FEM. 


105 


Figure  6-17  shows  a  screen  shot  of  this  simulation.  The  top  two  images  show  the 
temperature  distribution  at  the  front  and  rear  of  the  gage.  Although  these  distributions 
are  similar,  there  are  differences.  These  temperature  gradients  through  the  thickness  of 
the  beam  are  examined  in  the  bottom  left  panel.  At  this  measurement  point  in  time,  the 
variation  is  relatively  small,  less  than  0.1  K  at  any  point.  However,  a  time  history  of  the 
maximum  gradient  through  the  thickness  is  shown  on  the  bottom  right.  This  image  show 
that  as  the  shock  passed  over  the  gage  surface  and  for  a  small  period  of  time  afterward, 
the  thermal  gradient  reached  about  0.3  K  at  a  maximum. 
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Figure  6-17. 3-D  FEM  time  and  space  varying  heat  flux. 


This  simulation  proves  that  spatially  varying  heat  flux  could  have  contributed  to 
the  thermal  sensitivity  noticed  in  the  VT  supersonic  wall  shear  test.  This  combined  with 
the  possible  instrumentation  error  would 
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7  Conclusions  and  Recommendations 

This  research  provided  several  contributions  in  the  area  of  the  design  of  direct 
measure  skin  friction  gages  for  turbulent  compressible  flows  where  shocks  might  be 
present.  These  contributions  can  be  listed  as: 

1 .  Designed,  built,  and  tested  a  gage  that  measures  wall  shear  and  moment  (CF2) 

2.  Determined  most  suitable  gage  characterization  techniques 

3.  Developed  design  optimization  procedure  for  bellows-oil  configuration  gages 

4.  Designed  gages  for  short-duration,  high-enthalpy  testing  (CF3  and  CF4)  and 

examined  effect  on  design  space  for  blow-down  facility  and  flight  tests 

5.  Wrote  multiple  simulations  for  static,  thermal,  and  dynamic  analysis 

The  first  objective  was  to  build  a  gage  to  separately  measure  wall  shear  and 
moment.  CF2  accomplished  this.  In  a  pure-shear  environment,  where  voltage  output 
from  the  head  bridge  was  either  small  or  showed  signs  of  thermal  strain,  the  root  voltage 
output  provided  high-fidelity  results  for  wall  shear.  In  the  case  where  a  shock  was 
present,  CF2  captured  the  event.  Additionally,  the  oil-filled  configuration  completely 
eliminated  uniform  wall  pressure  as  a  source  of  apparent  strain. 

The  temperature  effect  seen  during  the  Virginia  Tech  Mach  4  pure  shear  case  was 
undesirable.  This  phenomenon  is  partially  attributed  to  the  long  transient  period 
experienced  for  tunnel  start.  In  fact,  it  is  possible  that  the  CF2  configuration  would  have 
worked  successfully  in  the  short  duration  testing  of  the  LENS  facility.  For  example,  the 
CF2  first  order  time  constant  was  about  0.34  s,  whereas  a  LENS  test  window  might  be  5 
ms.  Unfortunately,  the  CF2  gage  leaked.  This  combined  with  the  desire  to  attack  the 
thermal  issue  led  directly  to  a  complete  gage  redesign. 

The  need  to  develop  an  encompassing  design  methodology  first  required  the 
selection  of  tools  for  gage  characterization.  This  objective  was  met  through  the  careful 
examination  of  all  sources  of  strain  and  dynamic  response  characteristics.  As  a  test  of  the 
modeling  techniques  chosen,  each  of  the  four  gages  CF1-CF4  was  carefully 
characterized.  This  also  allowed  us  to  compare  the  gage  performance  to  the  nominal 
conditions  first  proposed  in  Table  1-4.  For  convenience,  the  expected  strain  output  is 
presented  again  in  Figure  7-1.  This  shows  that  despite  possible  sources  of  apparent 
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strain,  shear  output  dominates  the  signal  for  the  CF2  through  CF4  cases.  Uniform  wall 
pressure  is  not  even  a  factor  for  these  three  gages.  As  shown  during  gage 
characterization,  each  of  the  other  possible  sources  of  apparent  strain  provide  less  than 
10%  of  the  overall  signal  output.  Additionally,  all  of  the  gages  have  sufficient  dynamic 
response  to  settle  within  the  10  ms  steady  load  condition. 


Strain  Output 


□  Shear  ■  Inertia  □  Thermal  Gradient 
Figure  7-1.  Expected  strain  output  due  to  nominal  load  conditions. 

With  the  needed  modeling  tools  to  accurately  predict  gage  performance  in  hand, 
an  optimization  method  could  be  developed.  This  objective  was  successfully  met  by 
constraining  the  design  space,  virtually  building  and  testing  potential  gages,  and 
comparing  them  to  facility  constraints.  This  process  targeted  three  specific  test 
environments  and  considered  all  of  the  physics  expected  in  such  tests.  Direct  skin 
friction  gage  design  is  no  longer  a  matter  of  just  balancing  strain  output  with  vibration 
frequency  while  considering  inertial  and  temperature  effects.  It  is  now  a  full  optimization 
of  expected  physical  load  conditions  against  desired  output.  This  optimization  aspect  of 
the  research  may  prove  the  most  beneficial  for  future  work.  It  condenses  the  design 
process  to  a  set  of  state  and  design  variables  that  can  be  specified  for  a  particular  test 
environment.  This  combined  with  the  use  of  oil-filled  bellows  greatly  constrains  the 
design  process  and  allows  for  a  relatively  quick  turn  around  from  test  condition 
specification  to  gage  design. 
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The  impact  of  facility  specification  on  design  is  clearly  apparent  in  Table  7-1. 
This  table  shows  the  design  space  using  INVAR  for  three  very  different  test  conditions. 
The  high-enthalpy,  short-duration  facility  is  the  most  constrained,  ignoring  heat  flux. 
Additionally,  we  see  that  for  the  blow  down  facility  we  can  actually  consider  longer 
beam  designs. 


Table  7-1.  Optimized  design  space  for  INVAR  gages. 


Facility 


Design 

Space 


High  Enthalpy 
Short  Duration 


Blow  Down 


Flight  Test 


L 


>  th 


From  our  optimization  research  we  came  up  with  four  drivers  for  the  short 
duration  test  condition:  heat  flux,  minimum  frequency,  minimum  shear  output,  and 
maximum  acceleration.  But  let’s  look  at  how  test  time  affects  those  variables.  In  a  short 
duration  test,  we  need  quick  settling  time  and  even  high  heat  flux  has  little  effect. 
Conversely,  in  a  longer  test,  even  small  heat  fluxes  tend  to  provide  thermal  strain  issues. 
Figure  7-2  takes  a  second  rough  look  at  what  happens  to  design  space  for  two  very 
dissimilar  test  conditions. 
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Figure  7-2.  Design  drivers  and  effective  design  space. 


Our  objective  to  design  gages  for  application  in  LENS  was  met  through  the 
optimization  process  outlined  in  Chapter  5.  Additionally,  we  have  seen  the  expected 
output  from  these  gages  through  Gage  Characterization.  Although  these  gages  are 
currently  in  fabrication,  the  virtual  testing  by  computation  methods  has  shown  they 
should  meet  all  of  the  constraints  expected  at  LENS. 

The  final  objective  was  met  through  the  multiple  simulations  presented  in 
Chapter  6.  This  thorough  investigation  looked  at  dynamic  response,  thermal  response, 
and  instrumentation  effects.  Through  these  investigations  we  commented  on  possible 
sources  of  error.  Additionally,  we  were  able  to  perform  virtual  tests  and  look  at  expected 
gages  response  and  instrumentation  output. 

7.1  Short-term  recommendations 

The  design  optimization  section  provided  an  INVAR  and  aluminum  solution  to 
the  expected  load  conditions  for  a  specific  LENS  test.  These  gages,  CF3  and  CF4 
respectively,  are  currently  being  fabricated.  In  total,  4  INVAR  CF3  and  6  aluminum  CF4 
gages  will  be  built.  This  allows  for  several  different  strain  gage  configurations  for  each. 
The  quantity  and  type  of  configuration  being  built  is  shown  in  Table  7-2. 
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Table  7-2.  Gages  currently  being  fabricated. 


Gage  Type 

Quantity 

Strain  Gage 

Bridge  Type 

CF3 

2 

EA-00-03  ICE-350 

2  Half 

CF3 

1 

EA-00-15LA-120 

2  Full 

CF3 

1 

E  A-00-05 0- AH-3  5  0 

2  Full 

CF4 

2 

EA-13-03  ICE-350 

2  Half 

CF4 

2 

EA- 13-01 5LA- 1 20 

2  Full 

CF4  3-D 

2 

TBD 

4  TBD 

The  different  configurations  are  being  built  to  test  the  thermal  compensation 
proposed  in  the  theoretical  development.  The  half-bridge  configurations  work  like  the 
previous  CF2  gage.  The  full-bridge  configurations  are  intended  to  provide  thermal 
compensation  for  an  axially  varying  heat  flux.  Because  the  smallest  available  strain 
gages  have  lower  resistance,  a  slightly  larger  gage  is  also  being  used  in  a  full  bridge 
configuration.  This  allows  the  comparison  between  the  120  Q  and  the  350  Q  full 
bridges. 

Though  these  gages  have  been  designed  to  work  in  the  LENS  facility,  validation 
testing  will  be  conducted  in  the  Virginia  Tech  supersonic  tunnel.  The  test  results  using 
several  configurations  should  verify  if  an  axially  varying  heat  flux  is  a  source  of  apparent 
strain.  The  tests  should  also  show  if  a  full-bridge  configuration  can  provide  thermal 
compensation.  Additionally,  the  testing  should  show  whether  a  smaller  resistance  strain 
gage  can  provide  the  desired  output. 

The  expectation  is  that  the  INVAR  CF3  gages  will  perform  better  in  the  Virginia 
Tech  facility  due  to  thermal  characteristics  and  low  inertia.  Depending  on  the 
accelerations  in  the  LENS  facility,  the  CF4  gages  may  actually  perform  better  due  to  the 
lower  weight.  Also  depending  on  the  best  performing  CF4  configuration,  two  additional 
fully  3-D  gages  can  be  built. 

7.2  Mid-term  recommendations 

The  CF4  aluminum  gage  design  provides  a  unique  opportunity  for  additional 
strain  gage  configuration.  Because  the  beam  is  of  square  cross  section,  a  full  3-D  shear 
measurement  can  be  made.  This  would  merely  require  the  application  of  the  desired 
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strain  gage  configuration  on  the  opposing  beam  face.  This  is  the  intention  of  the  last  set 
of  CF4  gages  presented  in  Table  7-2. 

7.3  Long-term  recommendations 

There  are  myriad  approaches  to  attempting  skin  friction  measurements.  The  long¬ 
term  recommendations  in  this  section  provide  incremental  approaches  using  the  gage- 
type  developed  through  this  research.  The  cantilever  beam  type  direct  measure  gage 
cannot  completely  eliminate  the  effects  of  thermal  gradient  or  inertia,  at  least  not  using 
strain  gages.  It  is  possible  that  with  fiber-optic  technology  temperature  effects  could  be 
eliminated.  The  inertial  effects  though  small  can  also  be  greatly  reduced  for  a  cantilever 
beam  type. 

The  thermal  issue  is  important  to  consider  since  it  cannot  be  avoided.  This  is  an 
especially  important  consideration  for  flight  testing  where  heat  flux  acts  over  long 
periods  of  time  prior  to  measurement.  Research  using  fiber  optic  measure  of 
displacement  might  provide  a  suitable  alternative  to  strain  gages. 

The  evolution  in  size  from  the  CF1  through  the  CF4  gages  was  driven  partially  by 
choosing  the  FC-8-L  bellows.  It  is  possible  to  return  to  a  smaller  bellows  like  the  FC-3  in 
a  new  lightweight  design.  Though  the  assembly  process  would  be  more  complicated,  the 
reduction  in  weight  would  be  significant. 
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Appendix  A  Drawings 


Drawing  1.  CF3  element. 
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Drawing  2.  CF3  base  and  housing. 


116 


Drawing  3.  CF3  assembled. 
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Drawing  5.  CF4  base  and  housing 


119 


wmzm 


Drawing  6.  CF4  assembled. 
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